Unchained polygons and the N-body problem 
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Abstract 



We study both theoretically and numerically the Lyapunov families 

which bifurcate in the vertical direction from a horizontal relative equilib- 

I ) ■ rium in R'^. As explained in |CF1) . very symmetric relative equilibria thus 

give rise to some recently studied classes of periodic solutions. We discuss 
the possibility of continuing these families globally as action minimizers 
2 ' in a rotating frame where they become periodic solutions with particular 

symmetries. A first step is to give estimates on intervals of the frame 
rotation frequency over which the relative equilibrium is the sole abso- 
lute action minimizer: this is done by generalizing to an arbitrary relative 
equilibrium the method used in [BTj by V. Batutello and S. Terracini. 
In the second part, we focus on the relative equilibrium of the equal- 
Xj ' mass regular N-gon. The proof of the local existence of the vertical Lya- 

^^ , punov families relies on the fact that the restriction to the corresponding 

f*^ ' directions of the quadratic part of the energy is positive definite. We com- 

pute the symmetry groups Gr {N, k, rf) of the vertical Lyapunov families 
observed in appropriate rotating frames, and use them for continuing the 
families globally. 
^^ \ The paradigmatic examples are the "Eight" families for an odd num- 

ber of bodies and the "Hip-Hop" families for an even number. The first 
ones generalize Marchal's P12 family for 3 bodies, which starts with the 
k> ; equilateral triangle and ends with the Eight [CMi iMa2l ICFMI [CFll [S]; 

» I ■ the second ones generalize the Hip-Hop family for 4 bodies, which starts 

d \ from the square and ends with the Hip-Hop [UVl [CFll [TV] . 

We argue that it is precisely for these two families that global min- 
imization may be used. In the other cases, obstructions to the method 
come from isomorphisms between the symmetries of different families; this 
is the case for the so-called "chain" choreographies (see [5]), where only 
a local minimization property is true (except for A'' = 3). Another inter- 
esting feature of these chains is the deciding role played by the parity, in 
particular through the value of the angular momentum. For the Lyapunov 
families bifurcating from the regular A-gon whith A < 6 we check in an 
appendix that locally the torsion is not zero, which justifies taking the 
rotation of the frame as a parameter. 
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RELATIVE EQUILIBRIA 

AND THEIR LYAPUNOV FAMILIES 

1 Lyapunov families bifurcating normally from 
relative equilibria 

We consider relative equilibria rotating in the horizontal plane M^ x {0} C M'', 
that is, of the form 

where C is a central configuration, wi > is the frequency and J is the 
horizontal-rotation operator acting diagonally on the horizontal components by 
a rotation of 7r/2 and trivially on the vertical ones. Along the paper, various 
normalizations of C and various rotating frames will be considered, and the 
corresponding relative equilibrium will often be denoted by x{t). 

1.1 The horizontal and vertical variational equations 

The variational equations associated with a solution Xi = Xi (t) of the Newton's 
equations for N bodies in the Euclidean space (JS? , \\ ■ \\) 

mj- TTT, 2 = 1,...,7V 

are 

Along a planar solution Xi — Xi (i) (supposed to be horizontal) , the Pythago- 
ras theorem implies these equations split into horizontal and vertical parts. 
Namely, if 5xi = hi + Zi G C ® R is the decomposition of the variations in 
respectively horizontal and vertical components, (VE) is equivalent to the fol- 
lowing pair of equations: the complicated horizontal variational equation 



j-^j II J j-^j II J 



and the much simpler vertical variational equation 



i:™^^^' <"'^' 



i#i 



1.2 Vertical variations of a relative equilibrium 



As the mutual distances 



stay constant along a relative equi- 



librium motion, the corresponding vertical variational equation has constant 
coefficients: 
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OY z = Wz, where W is, up to a factor -2, the transposed of the Wintner-Conley 
matrix (i.e., up to a factor -2 and a transposition, it represents the endomor- 
phism A of the space of codispositions in [AC| ). The sum of the elements of any 
line of W is equal to 0. This implies that it acts on the space of dispositions 
T> = R^/(l, 1, . . . , 1)]R, to which z rightly belongs. 

The matrix W is symmetric for the mass scalar product, which means that 



z' ■ Y\Jz" = Wz' • z" 



where 



N 

^m,z[z'l. 



Hence the eigenvalues of W are real; because the Newton force is attractive, they 
are also negative (see |Mol Proposition 1] or [AC| V We will call their distinct 
values the vertical frequencies and denote them 



-c^?,- 



-w|,. 



-UJt 



e<N-i, 

where the LOk are chosen to be positive; note that by making W act on V we 
have taken away the eigenvalue 0. 

Now, let Zi, . . . , Z(7v-i) be a basis of V (which can be chosen orthogonal) 
consisting of eigenvectors of W with eigenvalues — w(l)^, . . . , —uj{N — 1)^, not 
necessarily distinct. The general solution Z{t) of (VVE) is of the form 

JV-l 

Z{t) = ^ Re iajZjc''^^^'''), aj € C, 



that is 



fe=i 
where each Wk is a complex eigenvector of W with eigenvalue —uj'^ 



1.3 What is known about the vertical frequencies 

First, one of the frequencies is that of the relative equilibrium, uji: it corresponds 
indeed to infinitesimal rotations around a horizontal axis. 

Let us now compare 071 to the other frequencies. Let /(x) = \x\'^ — "^miWxiW^ 
be the moment of inertia (i.e. the square norm in the mass metric), and let 
U{x) = X]j<7 \\x'-x\\ ^^ ^^^ potential function. Since a central configuration 



C is a critical point of the scaled potential 11 = 1^11, it is natural to write U in 
terms of U: 

dU{x)5x = -{x ■ 5x)I{x)-iu{x) + I{x)-^dU{x)5x 

and 

(fU{x){5x, 5x) = 2,{x ■ 5xfl{x)-iu{x) - 2{x ■ 5x)I{x)-^dtj{x)5x 

-\5x\^I{xyiu{x) + I{x)-^d^U{x){5x, 5x). 

Let now a; = C, the central configuration. Let us split a tangent vector 5x 
as before, into horizontal and vertical components h E C and z G R. Since 
dU{C) = 0, 

d'U{C)ih + z,h + z) = 3(C-/i)2/(C)~5t7(C)-(|/i|2 + |zn/-it/(C) 

+I{C)-id^U{C)(h + z,h + z). 

In particular, for a vertical variation, 

d^U{C)iz, z) = ~\z\^I{C)-W{C) + I{C)-^d^U{C)iz, z). 

But, taking the scalar product with x{t) = eJ'^'^^C of the identity 

VC/(C) = l(t) = -ujlx{t), 

where the gradient is relative to the mass metric, one gets I{C)^'^U{C) ~ ujf. 
Finally, one deduces that, if }\'{C)Zk = —uj\Zk, 

d^U{C){Zk,Zk) - {ujl~ijl)I-^\Zk\^. 

If we write any vertical variation Z = X)i=i ^i^(i) i^ terms of the orthogonal 
basis introduced in section [TT2l we get 

N-l 

r .2 _/_-\2\ 7-i_ 2l r. 

Hi) 



d'U{C){Z, z)=J2 K - ^i^f) i^^\Zi.- " 



Is was proved by Pacella (in the equal-mass case) and Moeckel (in the gen- 
eral case) that a planar central configuration C of at least 4 bodies is never 
a local minimum of U: more precisely, there exists always some Z such that 
d'^(J{C){Z, Z) < 0. This imphes the 

Lemma 1 (Pacella, Moeckel) For any central configuration C of at least 4 
bodies, at least one of the normal frequencies Wk is strictly greater than uji. 

This has a direct consequence on the Hessian of the Lagrangian action: let 

C(i) = Zfc(— t)=Re(W^fee'"^*). 



Then, for the action durmg time T = — , one has 

d^Aixit)mt),m)^n\Wk\'{L^j~L^l). 

whieh has the same sign as (PU(C){Zk, Zk). The proof of this formula is a direct 
computation: on the one hand, 

T 

on the other hand, we have seen that, for any t, dPU {C){C,{t) ^ C,{t)) = —ujl({t). 
This imphes the foUowing identity and hence the desired formula: 

d^u{x{t)){at),m)dt = -cjI r \at)\'dt = -o^iAWk?. 

Jo 

It follows that, if A^ > 4, a relative equilibrium is never a local minimizer of 
the action (compare |Clj ). 

1.4 Lyapunov families and their lifts 

_x(t),lift of 5 

Lyapunov orbit S 

/S0{2) 





Equilibrium E 
Relative equilibrium x{t), lift of E 
Figure 1 : Lyapunov family of a relative equilibrium 

The system may be reduced by fixing the angular momentum and quotient- 
ing by rotations around its axis; if the relative equilibrium x{t) is horizontal, the 
angular momentum vector will be vertical. Besides, by definition of a relative 
equilibrium, the solution x{t) projects to an equilibrium point E of the reduced 
system. Each pair of conjugate, purely imaginary eigenvalues of the lineariza- 
tion of the reduced vector field corresponds to a 27r/a;/t-periodic eigenmode. In 
good cases, these eigenmodes give rise to a family of solutions of the reduced 
system of period close to 2t: /cuk in the neighborhood of the equilibrium, called 
a Lyapunov family. 

In the simplest case where no other frequency of the equilibrium is an integer 
multiple of u (this is for example the case of the Hip-Hop family of the four- 
body problem), the Lyapunov center theorem [Mosj proves the local existence, 
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uniqueness and regularity of this family. In resonant cases, one can use either 
an argument involving a higher order normal form (such as in [CF2j for the 
Pi2 family), or a topological, existence argument (such as the Weinstein-Moser 
theorem |Mos| ) : for a detailed discussion on the simplest examples, see section^ 
Global existence theorems also apply [AYl ICMPj , but they give only general 
information on the nature of the families. 

Each Lyapunov family L lifts to the original phase space as a local one-para- 
meter (non-reduced, Lyapunov) family of invariant two-tori foliated by quasi- 
periodic solutions; on a given torus, any two solutions differ only by a rotation 
around the vertical axis (figure [T]). A non-reduced Lyapunov family associated 
with an eigenmode z{t) = KeZk e"^''' is tangent to the lift to the phase space 
of the cylindrical family of curves 

t ^ {x{t),Az{t~ If)), 

parametrized by a phase (p and an amplitude A. When observed in a rotating 
frame which puts into resonance the horizontal and vertical frequencies, these 
curves become a cylindrical family of loops, whose period is an integer multiple 
of 27r/wfe, say s2TT/u}k, s G Z \ 0. 

In a continuous family of rotating frames whose rotation frequency now 
depends on S (in general the rotation of the frame will vary, due to torsion, while 
it is fixed along the one-parameter family of solutions in the tangent cylinder), 
the quasiperiodic Lyapunov solutions x{t) themselves become periodic solutions 
whose homology class in T^ is indcpcndant of S. More precisely, if fJ,{S) £ R/Z 
is their monodromy, defined by 

x{t + 2t:/lu{S)) = eJ2'^^(^)/'^(^)a;(t), 

where uj{S) is the frequency of S and J, as before, is the rotation operator!^ the 

condition is that 

u(S) — zu ^ 

uj{S) 

The scaling invariance of Newton's equations (if x(t) is a solution, so is 
X~3x{Xt) for any A > 0) entails that each family can be freely rescaled. Fixing 
the norm of the angular momentum is a way to choose a normalization. But 
when moving away from the relative equilibrium, this may artificially lead to 
singularities, at elements of the family having zero angular momentum. Fixing 
instead the period (in addition to the direction of the angular momentum) is 
then a better choice of normalization, allowing to bypass zero-angular momen- 
tum elements. (On the other hand there is no hope with our techniques to 
bypass singularities where the period itself tends to infinity.) This allows to 
assume that in the above family of rotating frames the Lyapunov solutions all 
have the same period s2Tr/ujk- 



^Vertical colinear motions should be avoided since they correspond to a singularity of the 
reduction. 



When the loops in the cyhndrical family have some symmetries, the ques- 
tion arises whether the solutions in the the Lyapunov family, when observed 
in the above family of rotating frames, share the same symmetry. The unique- 
ness statement in the Lyapunov theorem implies that the answer is positive 
locally when the eigenmode is simple. However, the first eigenmode is always 
degenerate, which requires more analysis (see section H]). 

The really interesting question is to continue the family globally: if the 
integrated torsion effect allows continuation up to the inertial frame, we shall 
have completed the primary purpose of our study -to prove the existence of 
symmetric solutions in the inertial frame. If the symmetry group is rich enough, 
this program can be achieved using minimization of the Lagrangian action. This 
will be illustrated in section [5] with choreographic and Hip-Hop solutions arising 
from the regular A^-gon relative equilibrium. 

2 Minimizing properties of relative equilibria 

Assume here that the central configuration C is normalized so that e^'^C is a 
relative equilibrium. 

We study the action minimizing properties of the rescaled relative equilib- 
rium 

where 

C^ = (c. + n7)-ic. 

In a frame rotating with frequency zu, it becomes the T = 27r/a;-pcriodic loop 

Let A be the space of T-periodic loops in the configuration space with H^ 
regularity. For each value of ru we define the action Ar^{y) of y{t) G A as the 
action of the path x{t) = e^^*y{t) in the inertial frame: 

A^{y) = A{x) ^ J^ (^\x\^ + U{x)\dt^J f^\y + Iwy\^ + U{y)\dt. 
We call Au^ the action in a frame rotating with frequency va. 

2.1 Local minimizing properties 

In a frame whose rotation frequency is xu, a;^(t) has frequency lu. The La- 
grangian action during a given time T of a curve x{t) = {xi{t), . . . , Xn{t)) in the 
configuration space is 



A{x) = 



^\xit)\^ + Uixit)) 



dt. 



where, as before, |ip = ^mi||xi|p and U{x) ~ X]i<i \\T'^x \\ ^^ ^^^ potential 
function. As the kinetic and potential energies of x'^ are independant of time, 
its action during time T is readily computed to be 

A{C,u;,T,m) := A{x^) = {uj + m)i ^a, 

where a is the action of the normalized solution a;J(i) = e^*C. 

As expected, the action tends to zero when zu tends to —w. at the limit, the 
rotation of the frame alone accounts for the motion and the bodies have only 
to stay still at infinity in the inertial frame. 

We are interested in the local action minimizing properties under appropriate 
symmetry constraints of the members of this family in the following case: 



r 




27r 


27r 


-^k, 


T = 


= s — = 


- r — 


s 




L^k 


UJ 



where wi = 1, W2, • • • , Wfc, • ' ' ^I'c the frequencies of the vertical variational equa- 
tion associated with a;J(i) and r, s are mutually prime integers. Note that T is 
the minimal period of {x^{t), Zk{t)) in the rotating frame with frequency w, 
where Zk{t) has frequency ujk- 

Let us first compute the Hessian 

d'A{x){ts,) = / \\m^ +su{x{t)){m.m)\ dt. 

Jo '- ^ 

of the action when x = x^ and ^ ~ (0,0, Zk) with Zk a solution of (VVE) of 
the form Zk{t) — Re (VFfcC*'^''*) and Wk a complex eigenvector of WiC) with 
eigenvalue —ujI (see section 1.2). As J^ cos'^{uJkt)dt = J^ sin^ {LJkt)dt, we get 
immediately that 

\Zk{t)\^dt = L0l f \Zk{t)\^dt. 

Jo 



Hence 

d^A{x'^^){Zk,Zk) = I [ujl\Zk(t)\^+d^U{C^)iZk{t),Zk(t))] dt. 
Jo 

Now dU{x)X = X ■ Wx, where x and X belong to P (g) R^ = P^ ^nd W acts 
on each component. Since the mutual distances, hence W, stay constant under 
a vertical variation, we get 

d'UiC^)iZkit),Zkit)) - Z,it) ■ W{C^)Zk{t) - ~{C^/C)-^u:l\Zk{t)\\ 

that is 

d^A{xl){Zk, Zk) = {1-{lo + vjf)ul f \Zk{t)fdt. 

Jo 

In particular, d'^A{x^){Zk, Zk) = iff tu = ±1 — uj = ±1 — -ojk- This is easy 
to explain: A is the action in a frame which rotates in such a way that the 
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Not a local minimum 
Local miiiimimi" 




dMr«)(Zt, Zfr) < > V > " < " 



—uj — 1 —u) cj + 1 

Figure 2: Action of a relative equilibrium in a frame with rotation speed tn 

period of the relative equilibrium becomes s times the one of Zk (t) . It follows 
that Zk (t) defines a periodic Jacobi field and this implies that the corresponding 
variation belongs to the kernel of the Hessian of A. 

In the rotating frame, the solution x'^ possesses all the symmetries of the 
solutions of the (VVE) considered above. If the symmetry constraints that we 
impose are strong enough to allow only one type of solutions of (VVE) and no 
solution of (HVE), x'^ will be the only local minimizer of the action under thee 
constraints in the whole interval — (1 + oj) < ra7 < 1 — w. 

In the case of the regular n-gon configuration of n equal masses, this will be 
the case for 3 or 4 bodies; but already for 5 bodies, we shall see an example where 
the symmetry constraint does not discriminate between two different solutions 
of (VVE). In the next section we study the global minimization problem. 

2.2 Global minimizing properties 

This section, which as well as 15.11 is taken from the unpublished manuscript 
|C2| , develops an idea first used by BarutcUo and Terracini in |BT| . 



We shall now systematically use the shorter notation x{t) = gJ('^+™)*2; {^i 
place ofx'^it) = eJ'('^+^)*C^. 

Let G be a finite subgroup of 0(2) x S^v x 0(3). It has a natural action 
on the space A of T-periodic loops in the configuration space of the problem: 
g = (t, a,p) £ G acts on y{t) = {^yi{t),- ■ ■ , yAr(i)) £ A according to the diagram 



gy- 

The transformed loop is 



y: R/TZ x {!,..., TV} 
X {1,...,7V} 



Pi 



this convention is chosen so as to have a left action, i.e. {gg')q = gig'q)- 
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We suppose that y{t) = e~^'^*x(t) — e^'^^x belongs to the subset A*^ £ K of 
loops which are invariant under this action and we look for conditions on G and 
w which insure that the sole absolute minimizer of A^ is y. 

2.2.1 Choice of a new functional 

Following [BTj we are looking for an action functional y i— > ^^(y) with the 
following properties: 

1) Ar^iv) > A^iy) for any y{t) e A'=, 

2) Aixi attains its minimum value at y(t) 

3) If y{t) is a solution in the rotating frame, then Au,{y) — -^ro(y) if and 
only if y(i) = y(t). 

If such a functional exists, the sole absolute minimum of .4ro(y) is y{t). 
As in |BT| , from which the notations arc borrowed (but with different nor- 
malization), one looks for a functional of the form 



A^iy) = A{x) = 2 51 ^y^S + "^ I 51 ^''^^i. 



where, ii x{t) = (xi(t),--- ,XNit)) = eJ'^*2/(t), 

C^—/ \\x.{t)-x,{t)\\'dt= f \\y,{t)^y,{t)\\'dt 
Jo Jo 

and all the coefficients A, c, fiij are positive. 

2.2.2 Minoration of the potential part of the action 

By Jensen inequality applied to the convex function /(s) — s~^ , 

with equality if and only if ||xi(t) — a;j(i)|| is independent of t. 

Hence, for any choice of positive /iy 's, there exists a maximal c > such 
that 

r u{x)dt >T^J2 "^^^^ i^^^y^ >^tUY1 ^'^l^^o 

Indeed, as the function 
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is homogeneous of degree 0, it is enough to look for a constrained minimum of 
Si<7 f^if^j i^ij) ^ on X!«<7 f^ij^ij — I7 ^■^- there must exist a such that 

_2 

Vi < J, ^m,mj {(-jy^ = afi,,, i.e. V* < j, C" = l^^^j ■ 

Choice of the /i^: On chooses /i^ ~ fiij such that the f^ = I ^^^'^ 1 be 

equal to the squared mutual distances f? of the central configuration x (with 
relative equilibrium motion of period T in the rotating frame). In other words, 
one chooses 

p.ij = mimjf~y 

One checks that c = U^ , where U — U{x). We have proved the 

Lemma 2 



/ U{x{t))dt>{UT)i (^/2.,e 



with equality if and only if there exists /? > and S(t) G 0(3) such that for all 
t, x{t) = (3S{t)x. 

The only if part is because the mutual distances of the configuration x{t) 
being independent of t and proportional to the corresponding ones of S, this 
implies the existence of (3 and S{t). 

2.2.3 Minoration of the kinetic part of the action 

Integrating |x(i)p by parts between and T does not lead to a boundary term 
because, \i x{t) itself is not T-periodic, < x{t),x{t) > is T-pcriodic; one gets 

I \x\'^dt— I N^ mi||xi|pdt = — / y ^nij < Xi,Xi > dt. 
Jo Jo ,^1 Jo ,^1 

Let X = 2? (8) M^ be the quotient of (M^)^ by the action of translations, and 
A : A" ^ A" be defined by 

(Ax)i = ^ ^^(xj -Xj) = ^ ^(xj -Xj) if X = {xi,-- ■ ,xn). 
Then 



^m, / <{Ax)^{t),x,{t)>dt^ Yl rnimj_ 



"^ •Jj'i Jy 7 . •Jj'i ^ LLv 



Ju 7 *-*^'7, \AjL/ « 
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Definition. Let A^ be the smallest (eigenvalue) A such that there exists a 
solution y{t) G A*^ of the equation —y — 2Jvjy + vj'^y = XAy, or equivalently a 
solution x{t) = e^'^*y{t) of the equation —it = AAx with y{t) E A*^. 

Following BT], the Poincare inequality of the Kepler case is replaced by 
Lemma 3 For any y{t) E A^ , one has 

T 

with equality if and on if —it — A^Aa;. 

Proof. Define Aq as the minimum of the positive functional 



T 

\y{t)+iSwy{t)\^dt-- 


= f \iit)\'dt>\gy 

•^0 ^<J 



J{x) 



E.<,^T^/o'll^.W-^.WIPrf^ 



defined on the set of x{t) = e^^*y{t) in H^ with y{t) E A'^. The existence of 
a minimizer x{t) is insured as soon as the class of loops A implies coercivity. 
Writing the Euler equations, one gets that for such a minimizer, 

—X — J(a;)Ax. 



Lemma 4 Aq = A^. 

Proof. We have seen that if x{t) = e^'^^yit) with y{t) E A*^ minimizes J(x), 
it satisfies the equation —x = A Ax with A = J{x). The converse comes from 
the identity 



J(x) 



Y.i=l ^i /o < -*»' ^i(*) > dt 



Etl "^« Jo < {^X^){t),X,{t) > dt 

Being the minimum of a positive functional, A^ is necessarily non negative; 

this allows to write its definition as follows: 

Lemma 5 A^ is the smallest A > such that there exists a solution ^(t) of 
the equation — ^ = A^ such that y{t) = e"'''^*^(vAi) E A'~^ (and in particular is 
T -periodic). 

Remark If A is the space of all H^ loops of period T, Xmin — because of the 
lack of coercivity: the min is attained by motionless bodies at infinity. 

Interpretations of A: up to a harmless factor 2, a transposition changes 
A into the Wintner-Conley matrix associated with the central configuration x 
(cf. II. 2p . In other words, one checks that if AT = (A"i, • • • , ATjv) is a tangent 
vector at x to the configuration space (R'^)^ (or to its quotient by translations), 



dU{x)X ^ {X,-A3 
14 



Also, \i Z = {Zi, ■ ■ ■ , Zm) is a vertical variation, 

d?U{x){Z,Z) ^ {Z,-l^Z), 

which explains that A is symmetric with respect to the mass scalar product. 
Considered as living in the space M^ of vertical variations (resp. in its quotient 
by translations), the equation 

Z = -AZ 

is the vertical variational equation (VVE) at x{t). 

2.2.4 Minima of ^^^ 

Let g be the real function defined on M+ by 

The functional Aru is defined by 

\i<J 

The function g is strictly convex and its unique minimum is 

UT 



The minimum of Atu{y) = A{x) is attained when X]i<i P-ij^fj ~ ^min, that is 

i<j ^ij Jo (A^)3 

This implies the 

Lemma 6 A necessary and sufficient condition for the relative equilibrium y(t) 
to be a minimizer of the action Atu "in A^ , is that A^ = 1. 

Proof. If y{t) is a minimizer, J2i<j ™f™^ /o \\xi{t) — Xj{t)\\'^dt — UT a.nd the 

above identity becomes A^ — 1. 

Conversely, if A^ — 1, the unique minimum of the convex function g{s) = 
is + (C7T)5s-2 is attained at Smin — UT. But y{t) verifies Vi,j, ^fj = ^ijT- 
Hence J2i<j fiij^fj = ^^' which proves the assertion. 

Lemma 7 Let y{t) € A*^ be a solution in the rotating frame. Then, the equality 
•A.^{y) — -4ro(y) occurs if and only if A^ ~ 1 and y{t) = y{t). 
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Proof. -4ro(y) = A-!^{y) if and only if the equality holds for both the kinetic 
and the potential part of the action. Equality in the potential part implies that 
x{t) = e^'^*y{t) = (3S{t)x is a rigid motion. As it is a solution, it follows from 
[AC) that it is a relative equilibrium and hence of the form 

x{t) = (3Se'°''x, 

where 5' is a rotation. Now, equality in the kinetic part of the action implies 
that x{t) is a solution of the differential equation —x = XAx, where A = A^. 
Hence 

pa^Se'^^x = AA(/3S'e*"*x). 

But f3Se"^\xi, ■■■ ,xn)^ (/3S'e"*Xi, • • • , PSc^'^^xn) and hence 



Finally, 






Pa^Se^^'x = pSe'°"XAx = f3{Lo + myXSe'°"x, 

and 

a^ = {oj + w)'^\. 

So, x{t) = /3S'e'('^+^)^*x, and y{t) = e-'"*/?S'e*("+")^*i. But y{t) is 
supposed to belong to A'' and, in particular, it must be {T = — )-periodic, that 
is for all i. 

In particular, taking i = 0, this implies 

e'2-(^^^-f ) = 1, i.e. A = l, 

provided nj ^ —lu. This last condition is necessarily satisfied if y is a true 
solution, not lying still at infinity. From A = 1 and the fact that y and hence 
also a; is a solution, one deduces that /3 = 1 and hence y{t) — y{t). 

We conclude that 

Proposition 8 As long as X^^l, y{t) — xe*"'^* is the sole absolute minimizer 
of A„ in K^. 

Remarks 1) An elementary computation shows that 



-L / ^ 9 V^ -9 m . ,-~rm ^l 



: miTnir: 



■2 



= \0T, 
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while 



and 



minA^ = {Xg)H^ + \)UT = (AS)U^(y(t)). 

2) In the case of the choreography symmetry, that is of G = Z/iVZ acting by cir- 
cularly permuting N equal masses after one A^th of the period, this proposition 
is a rewording of BT]. It is proved there that it implies that, for any N, the 
regular A^-gon is the sole absolute minimizer of the action among choreographies. 

2.2.5 Not saying anything on the ItaUan symmetry 

As soon as relative equilibria with different central configurations are allowed 
by the group action, the method cannot work: a criterium which is essentially 
local cannot detect the difference in action between relative equilibria whose 
configurations are unrelated to each other. This is coherent with the fact that 
no conclusion can be drawn from the property that the criterion for absolute 
minimization does not apply. 

In order to illustrate this remark we look at the Italian symmetry 

x{t + -)^-xlt). 

It is the simplest invariance requirement which implies coercivity and it possesses 
two properties needed for our purpose: 

i) any relative equilibrium satisfies the Italian symmetry; 

n) any solution of the (WE) of a relative equilibrium solution satisfies the 
Italian symmetry. 

In order to compare easily with the formulae in the next sections, we consider 
loops in the configuration space which become 1-periodic in a frame rotating with 
frequency w. Moreover, when G — Z/2Z with the Italian action, we shall use 
the notation A^ instead of X^. 

Let x{t) — Q^(^'^+'^)tx be some relative equilibrium {x can be any central 
configuration). Let a;i,W2,''' be the vertical frequencies (i.e. the ones of the 
Vertical Variational Equation —z = Az) of a relative equilibrium with similar 
configuration and frequency oji (see [CFlllMo] ). 

Lemma 9 Whatever the central configuration x, and whatever the dimension 
(two or three) of the ambiant space, whatever the value of vo , 

V^<inf^. 
V fc Wfc 

Corollary 10 The criterion for being an absolute minimizer does not apply to 
the Italian symmetry as soon as N > A: A^ is always strictly smaller than 1. 
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Proof of the lemma. We shall denote hy lui — 2tt + zu, lj2,- ■ ■ the vertical 
frequencies associated with the relative equilibrium x{t) — cce*'^^'^^'^-'*, i.e. 

CJfe = — [Ztt + ZU). 

Let us consider a horizontal solution ^(i) of the equation — ^ = A^, of the 
form (we identify the horizontal plane to the complex plane) ^(t) = Wke^'^''t. 
The periodicity condition of y(t) = e*™*^(vAi) and the Italian symmetry are 
satisfied if and only if it is of the form y{t) = Wfee™^'^** with m odd, i.e. if 
'i'fevA — ZU = m2TT, that is ^^(27r + vj) — m27r, that is 

r— UJl I ni2TT + ZJ 



LUk \ 2tt + w 

One concludes by choosing m = I. 

Proof of the corollary. Let Z{t) = Re (VFjte*"'=*) be a solution of the 
(WE) equation associated with the relative equilibrium x{t) = Axe'"^* with 
configuration similar to x and frequency ui . Here, Wk is a complex eigenvector 
with eigenvalue uj^ of the associated A operator. Finally, let 

at)^Z{^t)^Re{Wke'^''). 

^k 

Then, for the action during time T — —, we have proved in 11.31 that 

d'A{x{t))iat),at))^^\Wk\'{cj!-u;l). 

The corollary follows from the Pacella-Moeckel lemma recalled in 11.31 

Remark 1 Not knowing the result of Pacella-Moeckel, one could have deduced 
that, in order that the relative equilibrium solution x{t) = xoe''"^* be the sole 
absolute minimizer of the action in A^* , it is enough that it be a local minimizer 
with respect to the vertical variations, which would have been a mirific result 
indeed if it had not been an empty one, except in the case of the equilateral 
relative equilibrium of 3 bodies. 

Remark 2 In the following sections the choice of the groups Gr/s {N, k, rf) will 
eliminate the non trivial horizontal variations, since the only horizontal loops 
which are invariant are relative equilibria of the regular iV-gon, possibly with 
a different frequency and a reordering of the bodies. The same was true in the 
case of choreographies, studied in [BTj : the regular A^-gon is the sole central 
configuration whose relative equilibrium motions are choreographies. 

The main difference with the study below of the groups Gr/siN,k,r]) lies in 
the fact that the two integers m and k will not he any more independant! 
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THE CASE OF THE REGULAR iV-GON 

CA UTION: It will he convenient from now on to number the bodies from to 
N —\ considered as elements ofL/nT,. Moreover, we suppose that the mass of 
each body is equal to 1. 

3 Infinitesimal continuation 

3.1 Vertical variations 



We focus on the central configuration C — (1,C,...,C), with A^ > 3 and 
^ = e*^'^/''^, in the horizontal plane identified with C, and on the associated 
relative equilibrium 

x{t) = (xo(t), Si(i), . . .,XN-i{t)) = e^'"!* C 

of the equal mass A^-body problem. 
Let 

Pk=r^^^+k = \C'' -II k=l,...,N-l. 

In particular, pN-k — Pk ^i k < [N/2]. The vertical variational equation (VVE) 
reads 

-^ — 



( ^0 \ 



f-^i 



Pn-1 



Pi 



\ Pi 



1 



P'i 
1 



1 



N -1 



P%-2 



( ^0 \ 

Zl 



-Ej.) v--^J 



Such a circulant matrice has an explicit basis of eigenvectors: a basis of complex 
eigenvectors is formed by the 



Xk = ic',e', 



("" = !), k = 0,...,N-l, 



with eigenvalues 



Pn-1 



(W-l)fe ^ _ 



E 



I _ pk 



i<j<N-i ^i 



Pj Pi P2 

In particular, Ai — —uj\ corresponds to the frequency of the relative equilibrium. 



The cases N odd and N even behave slightly differently from one another. 

(i) Case N ~ 2n+ 1 odd: the 4n-dimensional phase space T)^ of the varia- 
tional equation splits into the direct sum of 2n invariant cigenplanes. Indeed, 
there are exactly n distinct mutual distances 



Pi = PN-l, P2 = PN-2, 



, Pn — PN-n = Pn+1 
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(in increasing order: see lemma [TT] below) and C,^ — C}-'^ ''\ Hence, for k = 
l,...,n, 

(the eigenvalue Aq = 0, with eigenvector Xo = (1, 1, . . . , 1), disappears after the 
quotient by vertical translations) . 

To each eigenvalue Afe (1 < /c < n) of the circulant matrix, corresponds 
the 4-dimensional space of solutions of the (VVE) consisting in solutions of the 
following form, where the amplitude A and the phase to are parameters and 
^k = v— Afc : 

z{t) = ('ARe(e*'^''(*-*«)),ARe(C''e*'^''(*-*«)),...,ARe(C''^^"^^e'"''(*^*o))] 

and 

z{t) = (ARe(e*"''(*-*°)),^Re(C''e*"'=(*-*»)),...,^Re(C''(^-i)e*'^''(*-*«))) . 

(ii) Case N = 2n even: the 4n — 2-dimensional phase space splits into the 
sum of a 4 (n — l)-dimensional invariant space in which the solutions are of the 
same form as in the odd case and an invariant plane which corresponds to the 
eigenvalue A„ of the circulant matrix. The eigenvalue A„ plays a special role for 
the sole reason that C,^ — —1 is real and hence corresponds to a 2-dimensional 
space of solutions and not a 4-dimensional one. 

3.2 The vertical eigenvalues 

Partly following the method of [PW[ IMo| . we prove the following fact, which 
will be used in sections |3] and [5] 

Lemma 11 The half sequence {Xk)i<k<N/2 's negative and decreasing. 

Case of an odd number of bodies: iV = 2n -I- 1. Using the fact that 
(^~i = (i and pN-j = Pj, and introducing 9 = 2tt/N, we get 

1 v-^ 1 — cos jk6 

4 sin - — 

j=l °"^ 2 

Three successive discrete derivations yield 

6X, = \k+i-Xk = -\Y:U ^S^^^ (l<fc<n-l) 

6^\u = ^Afe+i-5A, = -E;=i^2LZi^ (i<fc<n-2) 

5^\k = S^Xk+i-6^Xk=2j:;^^smj{2k + 3)l (l<fc<n-3) 

_ ^ sin((2fc+3)(n-|-l)|)sin((2fc-|-3)»|) 
^ sin(2fe-|-3)f ' 
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We want to show that Xk decreases with k, or that SXk is negative. The 
reason for derivating three times is that S^Xk is a trigonometric polynomial and 
that its sign is related to the convexity of 6Xk ■ 

Notice that these sequences still make sense for larger values of k, provided 
that the denominators do not vanish. Besides, the property that 5Xk is < for 
extremal values of k is easier to check outside the initial interval of definition: 



1 1 



J=l 



(because jf = J 2n+i ^ "") ^^'^ 



^ ~ sm 7 7j 



smjTT 

2 



0. 



Hence it suffices to show that 6Xk is convex with respect to k, i.e. S^Xk > 0, 
over {0, ...,71 — 2} (see figure [31). 



u 


' 


' 


' 
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_:¥' 
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_..-H'' 


- 
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'"■+---■ 


----■■^"'' 







1 2 3 4 5=n 

Figure 3: Graph of 6Xk for n = 5 (2n +1 = 11 bodies) 

In the above expression of S^Xk, the denominator is the sine of 

,e 2n-l n TT 

0<(2fc + 3-<- TT < 77 if 0<fc<n-2, 

^ ^4-2n + 12 2 -- 

and thus is > 0. It remains to check that the numerator 



j/fc = sin ( (2fc + 3)(n + 1)- 1 sin ( (2fc + 3)n- 



itself is > 0, which follows from linearizing: 



^fc = t: cos — 

Z n 



fc + I TT 

2->0 if 0<k<n-2. 



i 2 

2 ^ 
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Case of an even number of bodies: N = 2n. In the sum defining A^, all 
terms come pairwise (indexes j and iV — j), except that of the n-th body. Let 
ej — 1/2 if j — n and ej = 1 otherwise. Then 



1 x-^ 1 — cos jkO 

-> fi 3-z — (l<k<n). 



4" sm' 2 



It follows that 

c-, 1 v^n sinj(2fe+l)f n ^ 1 ^ i\ 



sm'^ J 7 



\k = 



-E;=ie.™7f^ (l<fc<n-2) 

2E"=ieiSinj(2fc + 3)| (1 < fc < n - 3) 

_ ^ sin((2fc+3)(n+l)|)sin((2fc+3)n|) , ^ -, ^fc 
^ ^ sin(2fc+3)f ^ '^ ^'' ■ 

Again, we will check that SX^ is < at the boundary of some interval of integers 
containing {1, ..., n — 1}, and that SXk is convex inside this interval. 

Similarly to the odd case, obviously we have SXq < and 5Xn-i/2 = 0. 
Unfortunately, the latter equality is not of any direct use since the index is not 
an integer. Let us rather show that 




is < 0. It suffices to see that 



is < 0. Note that in these sums if j < n — 1 is odd the pair of terms of indices j 
and j + 1 have a negative contribution. Hence, if n is odd, the last but one given 
expression of a„ shows that a„ < —1/2; and, if n is even, the same estimate 
follows from the last expression. 
It remains to show that 



(5'Afc = 



2 sin {{2k + 3)(n + l)f ) sin {{2k + 3)n|) + (-1)'= sin(2fc + 3)^ 

k ^ 



sin(2fc + 3)| 

is positive ifl<fc<n — 2. The denominator being > 0, focus on the numerator 
i/fe. Setting a = {k + |)^, we get 

i^fe = 2 sin(n + l)Q:sinnQ! + (— 1) sina; 
splitting [n + l)a into na + a and partially linearizing yields 

( \ 

cos a (1 + cos2nQ:) + tanQ; (sin2na + (— 1)*^) > 0, 

which proves the result. 
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3.3 The symmetry group Gr/s{N,k,ri) 

We will now analyze the symmetries of vertical eigenmodes in the full phase 
space. Recall that most of the time a vertical eigenfrequency is degenerate 
and has vertical multiplicity 4. We will let the amplitude and the phase be 
respectively A = 1 and io = 0, and consider the most symmetric generators of 
its eigenspace i.e. 2-frequency motions of the form 

2^jW== (C^e*'^i*,Re(C''''^'e''^''-*)) e C x K, (j = 0, • • • , TV - 1), S{N,k,T]) 

where rj = ±1 and, like in the former section, ( — e*"" and k G {1, ...n}. When 
observed in a frame rotating around the vertical axis with angular velocity zu 

such that 

Wi — 07 r ^ 

— = - e Q, 

LUk S 

the horizontal and vertical frequencies are set into resonance and the motion 
becomes periodic of period T = 22L£; 

a;,(t) = (C^e^5-^*,Re(C'"^^"e^'^'=*)) (j - 0, ..., TV - 1) Sr/s{N,k,r^). 

The discrete symmetry group of such a motion is seeked as before as a 
subgroup of 

Go = 0(M/TZ) X Y.{N) x 0{R^), 

where g = (r, cr, p) G Go acts naturally on the space of T-periodic loops: 

Let Gr/siN,k,r]) be the stabilizer of Sr/s{N,k,ri). The group structure of 
Gr/siN,k,ri) does not depend on r and it depends on s only through the fact 
that we are looking to the solution during a time interval s times longer than 
the minimal period of the relative equilibrium in the inertial frame. We let (s, k) 
be the gcd of s and fc, and s = (s, fc)s', k = (s, k)k' . 

Lemma 12 Gj,/siN,k,ri) is a semi-direct product of an Abelian group H of 
order 2Ns by Z/2Z. The group H is an extension by Z/2Z of a group K which 
is itself an extension o/Z/iVs'Z by Z/(A:,s)Z. 

Proof. 

1° Restriction to a subgroup G I of Go- Certainly elements of Gr/s(-^, fc, ??) 
stabilize the regular A^-gon relative equilibrium (as the horizontal component of 
any infinitesimal vertical variations), as well as the cylinder of infinitesimal ver- 
tical variations. Hence Gj./s{N, k, rj) is contained in the subgroup Gi consisting 
of elements g = (r, a, p) € Go satisfying the following conditions: 

— The isometry p e O(M^) is of the form p = [phor, Pver) : C x M h^ C x M 
where 

Phor{h) = e''^''°'h or e^^Tra^ ^^j p^eriv) = C^'^W 

with a e M/Z and f3 e Z/2Z. 
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— If we set ^ = ±1 according to whether phorih) = e'^'^'^^h or e^'^^°'h, 
r-i(t) = ^{t - 9) with 9 e M/TZ. 

Hence an element g (z Gi can be identified with a quintuple 

(6',CT,a,/3,0 e R/TZ x Y.{N) x K/Z x Z/2Z x Fa, 
the multiplication law being given hy g — g' g with 

9 = 9' + ^'9, 5- = o-'cr, & = a' + ^'a, /? = /?' + /3, ^ = C'6 
Jn i/ie rest of this proof, let us set the time unit so that LUk = 27r, i.e. T = s: if 

X e Sr/s(N,k,l]), 

_ / exp (27^^ [a + j^^a-^j) + f (t - 9)]) \ 
5^^' - 1^ cos (27r [f + j^ka-^j) + ^(t - 0)] ) ) ' 

2 ° Restriction to a subgroup G2 ofGi. Setting S = (,(7^^{j) — j G Z/A^Z, 
the symmetry equation gx = x reduces to 

a + jf- ^9 = (mod 1) 
^^ + kr]^-9 = (modi). 

The first equation shows that 5 is independent of j, i.e. S^a is a circular permu- 
tation. Hence, using the fact that S^/3 = /3 (mod 1), 

c^^^O-i (mod 1) 

61 = I + kr]-^ (mod 1) 

These equations completely determine a G R/Z as a function of 5 and but 
G M/sZ is only determined mod 1, as a function of (6, /3, £_) E Zn x Z/2Z x F2. 
Let 

G2 = (R/sZ X Z/7VZ X Z/2Z) x F2 = {{9,6, P,^)}, 

where the semi-direct product is defined by the law g = g' g, with 

9 = 9' + C'O, S = 5' + ^6, p = f3' + (3, i = ^'^ 
Then Gr/s{N, k^rj) identifies with the subgroup of G2 defined by the equation 

= h kri — (mod 1). 

3 ° Group structure of Gr/s{N,k,rj). 

(i) The semi-direct product structure of G2 goes down to Gr/s{N, k, rf). Indeed, 
the section ^ 1— > (0,0,0,^) of the projection G2 -^ F2 takes its values in the 
subgroup Gr/s{N,k,r]). Hence Gr/s{N,k,r]) = H xi Z/2Z, where H is the 
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subgroup of the abelian group R/sZ x Z/A^Z x l^jTL defined as the set of triples 
(6*, (5, /3) which satisfy 

e=- + kri— (modi). 



('ii) Let K be the kernel of the group honiomorphism {9, 5, j3) ^-^ (3 from H to 
Z/2Z. It can be identified with the subgroup of R/sZ x TLjNTL defined as the 
set of pairs (0, b) which satisfy 

9 = krj— (mod 1). 



('iii) One has the exact sequence 

-^ Z/Ns'Z -^ K ^ Z/(fc, 3)1 -^ 0, 

where the first arrow sends 1 to {krj/N, 1) and the second one sends {9, 5) to the 
class oi 9 — krjS/N. This ends the proof of the lemma. 

The precise structure of G,./s(iV, fc, 77) depends on s, N, fc, 77. We study it in 
some cases. 

Lemma 13 // one of the following conditions is satisfied, H is isomorphic to 
K X ■L/2Z: 

1) s is odd 

2) s and N are even, k is odd. 

Proof. 

1) If s is odd, the mapping f3 1— > (^,0,/3) defines a section of the projec- 
tion from H to Z/2Z. Hence the mapping (9,6,0) i-^ ((0 — ^, (5),/3) is an 
isomorphism from H to K xl/TL. 

2) If s = 2(T is even, a section must send 1 onto an element of order 2 in H 
not belonging to K, i.e. an element of the form {9, 6, 1) such that there exist 
a, 6 G Z with 

2(5 = bN, 1 + 2fc?7— = 1 + 6/C77 = as = 2aa. 

From the second equation, it follows that b and k must be odd. The first one 
then implies that N is even. If all these conditions are satisfied, /3 1^ (/3|, /S^i P) 
defines a section and hence the mapping (9, 6, 13) ^^ ((^ — /3|, (5 — (3^), 0) is an 
isomorphism from H to K x'L/2'L. 

Lemma 14 // (fc, s') ^ 1, K is isomorphic to l/Ns'Z x Z/(fc, s)Z (and hence 
to Z/NsZ if (Ns', [k, s)) = I); if [Nk' , [k, s)) ^1, K is isomorphic to Z,/NsZ. 
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1) If (fc,s') = 1, the mapping which sends 1 € 'L/{k^s)'L to (s',0) G X is a 
section. 

2) If {Nk\ {k,s)) = 1, this imphes (N, (fc,s)) — 1, and hence the existence 
of integers I' and 6 such that A^l' + riS{k,s) = 1. Let us set / ~ I'k' and 
9 = krjS/N + I. We assert that the element {9, S) is of order Ns. Indeed, 
p{9, 6) = E K if and only if there exists integers a, b such that 

pkrjS/N + pi = krjb + pi = as, p6 ~ bN. 

From its definition d satisfies {N,6) — 1. Then, the second equation above 
implies the existence of an integer p' such that p — p'N, hence 6 = p'5 and the 
first equation above becomes p'{{k, s)riS + Nl')k' = as, that is p'k' = as. As 
{Nk' ,{k,s)) = 1 implies (fc',(fc,s)) = 1 and hence (fc',s) = 1, there exists an 
integer p" such that p' — p"s, and hence that p = p" Ns. 

Corollary 15 When s = 1, Gr/s{N,k,rj) is isomorphic to the direct product 
Dn X Z/2Z, where Dn, of order 2N , is the dihedral group. In particular, when 
N is odd, it is isomorphic to D2N ■ 

3.4 Invariant loops 

An s-periodic loop of configurations x{t) = [Xi{t), ■ ■ ■ , x^it)) is invariant under 
the action of Gr/s{N, k, rj) if and only if, for every {6, 5, f3, ^) e G2 representing 
an element of Gr/s(N, k, rj), i.e. such that 9 — '^ — krj-^ = / G Z, one has 

V? e Z/A^Z, x,{t) - px5(,+5) (e(t - 9)), 

where the action of p on M'^ = C x M is defined by 

p{h, z) ^ {e''^''°'U,e'''^z) with a = -9-— (mod 1), 

where /i« = /i if ^ == +1 and /i« = /i if ^ == -1. 

Before looking at remarkable classes of invariant loops, let us make some 
general comments: 

1. Taking (6* = i + /, 5 = 0, /? == 1,^ = 1) and setting Xj{t) = {hj{t), Zj{t)), 
we get 

h,{t) = e*2.;(i+0/,^.(t -^+1), z,(t) ^ -z,{t - i + 0- 

In particular, if s = 2/ + 1 and r are odd, an invariant loop possesses the Italian 
symmetry. This is the case of the Hip- Hops and of chains with an odd number 
of lobes (see below). 

2. Taking {9 ^ 0,5 = 0, (3 = 0,£, = -1), we get 

hj{t)^h-ji-t), z,{t) ^ z^,i-t). 

In particular, when t = 0, the configuration is always symmetric with respect 
to the vertical plane containing the first coordinate axis. 



26 



3. Taking (61 = ^, (5 = 1, /3 = 0, ^ = 1), we get 



h,(t) = e^^-^h,+i{t~^), zj{t) = z,+i(t- ^). 
where A = -krj — 1. 

More generally, the smallest non zero value of 6* = |+fc7/;|+/ = N(p+2^)+2k7iS 
is 6*0 = 2jv because P + 21 and 6 are arbitrary integers. 

One can distinguish two cases in the action of G^/siN, k, rj) on the space of 
loops of A^-body configurations. In the first one, exemplified by the P12 family 
for 3 bodies, no a priori spatial symmetry of the configuration exists for all 
times, and even a = for any element of the group. In this case, the group acts 
on the full space of similitude classes of A^-body configurations. The condition 
is that for any / e Z, (S e Z/NZ, /3 e Z, f [f + kr]-^ +l]-± =0 (mod 1). This 
is equivalent to s = 1, r = 2r', 2r'krj— 1 = (mod A^) (which implies that A'' is 
odd). The case - = 2r' = A^— l,fc = 1,77 = — 1 corresponds to the "unchained 
polygons" with A^ — 1 lobes. 

In the second one, exemplified by the Hip-Hop family for 4 bodies, such a 
symmetry does exist: there is a group element such that ^ Q, 5 = —1, a = 
1/A^. In this case, the group acts on a subspace of the space of similitude 
classes of A^-body configurations. The condition is that f ~ 7v ~ '^ (mod 1) 
which, because fc < y) implies /3 = 1 and k-q ^ ^ (the choice of 77 — ±1 is 
then immaterial because e^^ = e~'^^ , so we set 77 = 1). The A^ = 2A^'-body 
Hip- Hops studied in jTV[ IBCPSj fall into this category; the group element is 
{6 = 0,5 ^ —!,£, = 1,q: = jf), which sends one body onto the next one by a 
rotation of ^ followed by a change of sign in the vertical component. We now 
study some types of symmetric loops in the two cases. 

The choreographic symmetry (bodies gather by pairs, at least, following 
the same curve in space): this corresponds to a group element g with ^ — 1,6 ^ 
0,a — 0, l3 ~ 0. The equations defining the symmetry group become 

i-'-e^O (mod 1) 
kr]j^-e = (modi). 

Choreographies can be simple or multiple. They are sim,ple if all the bodies 
lie on the same curve, i.e. if in addition a has a unique cycle: 

{S,N)^1. 

Lemma 16 Sr/siN,k,ri) is a simple choreography if and only if 

s — krjr = (mod A^). 

Proof. We rewrite more explicitely the conditions for being a choreography: 
there exists integers l,m such that 

s5 - NrO = IsN 
kr]5 - N0 = mN. 
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This is equivalent to 

(s — k'qr)5 — (Is — mr)N 
k-qS - N9 ^ mN. 

Due to the condition of simple choreography (S, N) = 1, this is equivalent to 
the existence of an integer a such that 

Ls — mr = aS 
s — krjr = aN 
k-qd - N9 = mN. 

This proves the "only if" part of the lemma. 

For the "if" part, let us suppose that there exists an integer a such that 
s — krjr — aN. As (r, s) == 1, we can choose I, m such that Is + mr = a. We 
then define 

9= (e=^-m,S=l,^=l,a = 0,p = o) e Gr,s{N,k,ri), 
which completes the proof. 

Remark The condition of simple choreography could of course have been 
obtained directly from the formulae for invariant loops or by writing that for all 
j, one has (^^'''^^e*'^''*) " — ^^6* = "'=*, which expresses that all the bodies belong 
to the curve described by the first. 

Corollary 17 For any N,k,ri, the set of values ofr/s for which Sj,/s{N,k,r]) 
is a simple choreography is dense in K. 

Proof. The integers N,k,ri being given, r,s,a must be such that 

(r, s) = 1 and s — krjr = aN. 

The primality condition is equivalent to the existence of integers I, m such that 
Is + mr — 1 , that is (m + lkri)r + l{aN) — 1. As (m + lkri, I) can be an arbitrary 
pair of integers, this is equivalent to (r,aN) — 1. But the density of the r/s is 
equivalent to the density of s/r — krj + aN /r and hence to the density of the 
aN /r or equivalently of the r / aN , with the unique condition that (r, aN) — 1. 
As the irreducible fractions of the form r /N"^ are dense, the corollary follows. 

Eights and maximal chains symmetries The two choreographic cases with 
extreme values of k are G2{N = 2n + 1, n, —1) (which we call the Eight sym- 
metry) and Gn-iIN,!, —1) (resp. the maximal chain symmetry). In the first 
case, the group element (9 = jj^,d = I, (3 — 1, ^ — I) leads to 

Xj{t) =Xj+i{t~—), Zj{t) = -Zj+i(i- — ); 
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in the second case, the group element (6 = 2^, (5 = ?i, /3 = 1, ^ = 1) leads to 

Xj{t) = Xj+nit - — ), Zj{t) = -Zj+nit - —). 

Solutions invariant under such a symmetry are called unchained polygons be- 
cause the horizontal rotation of the regular A^-gon is unfolded in the vertical 
direction, into a choreographic chain with 2 and A^ — 1 lobes respectively. Ex- 
amples are shown on figures [H 2 and 21 1. 

The Hip-Hop symmetry Writing k — {N,k)k, N = {N,k)N, we notice 
that 

^vkji ^ ^vkj2 jf ^^^ o^iy jf j^ _ j^ jg ^ multiple of N. 

Hence the bodies are divided into A^ regular {N, fc)-gons which remain horizontal 
at each instant. Independantly of the chosen rotating frame, the motion is 
invariant under a group element g of the form 

If we make the simplest possible choice r/s = 1, we get a Hip-Hop motion. The 
name is a reference to the original Hip-Hop solution where the two diagonals of 
a square stay horizontal while undergoing vertical oscillations (see figured 6). 
Note that this symmetry can be combined with a choreographic symmetry. 

A few examples of Lyapunov cylinders are shown on figure [H They are re- 
spectively tangent at the relative equilibrium to the following Lyapunov families 
(the terminology will be explained in section [S]) : 

1. 5*4(5, 1, —1) Maximal unchained polygon {£ — 4) 

2. 5*2(5, 2, -1) Five-body Eight 

3. 5*3(4,1,-1) Unchained polygon (Gerver family) 

4. 53/2(4,2,1) Terracini-Venturelli choreography (Hip-Hop family) 

5. 5i(6, 1, -1) Six-body Hip-Hop 

6. 52(6, 2, 1) Yet another six-body Hip-Hop. 

3.5 Isomorphic symmetries 

In order to better understand minimizers of the Lagrangian action among loops 
which in the rotating frame are invariant under the Gr/s{N,k^rj) action de- 
scribed above, we study the cases when the actions of two groups Gr/s{N, k, rf) 
and Gr' /s'{N, k' , rj') coincide up to a relabelling j i-^ j' = I](j) of the bodies. A 
first condition is that the periods coincide, i.e. s' = s. 

We will stick for simplicity to the case s = s' = 1. 

Let as usual C = ^^■^•'■/N- the horizontal relative equilibrium whose j'th body's 
motion is ft.'/(i) — C"* e^'^*'' * is G'r'(A^, fc',77')-symmetric. Hence a necessary 
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Figure 4: Some examples of first order solutions having special symmetries in a 
resonant rotating frame. 



condition for the two actions to coincide up to the relabelling S is that the 
horizontal relative equilibrium whose jth body's motion is 

be G'r(Af, fc, ?7)-synimetric. This means that for any (5,(3,^) G Gr{N,k,ri), 

Using the identities 9 = § + krjj^ (mod 1) and a = -6 — j^ (mod 1), this 
becomes 

ir~r')[^ + kv^]-^ + ^m^ij+6))-nj)]^0 (modi). 
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5 — 0, (3 ^ 1,^ = 1, wc get a first necessary condition 

1) r -r' = 2pe 2Z. 

Taking then S = Q, (3 = 0, we get ^S(^j) = S(j) (mod N). In particular, if 
^ == — 1 and j = we get E(0) — (mod N). Finally, taking S = 1,(3 = 0, we 
get E(j + 1) - T,{j) = 1 - 2pfcr7 (mod N), hence I](j) = (1 - 2pkr])j (mod TV). 
Interchanging the roles of the two groups we see that necessarily S~^(j') = 
(1 + 2pk'ri')j' (mod A^). This yields a second necessary condition 

2) (1 - 2pkr])(l + 2pkri') = 1 (mod N). 

It turns out that these necessary conditions are essentially sufficient. Indeed, it 
is enough to replace the condition 2), which can be written — 2p[— fcry + k'rj' — 
2pkrik'ri'] = (mod A^) by the slightly stronger condition 2') below: 

Proposition 18 A necessary and sufficient condition for the existence of a 
permutation & of Zi/Nlj such that for any Gr'{N,k',ri')-invariant loop x'{t) 
whose j' th body motion is (/i',(t), z', (t)), the loop x(t) whose jth body motion is 
{hj{t), Zj(t)) := (hL,Jt),zL,Jt)) be Gr{N ,k,rj) -invariant is that the following 
conditions be satisfied: 

1) r - r' ^ 2p e 2Z, 

2') -kri + k'rj' - 2pkr]k'r]' = (mod N). 

The permutation & and its inverse are then respectively j i-^ j' ~ {l — 2pkri)j 
and f ^ j = {1 + 2pk'-q')j' m Z/NZ. 

Proof. We already know that the asserted form of the permutation 6 is nec- 
essary. What remains is a direct computation: in order to show the invariance 
of the loop x{t) under the element of Gr{N,k,r]) represented by (5, /3,^), one 
must find an element of Gr'{N,k',ri') repesented by {S',(3',£,') such that 

and 

Straightforward computations show that this reduces to condition 2'), and 
that the unique solution is 

(5' = (1 - 2pk7^)s, (3' = f3,e = e, 



Examples We note G = G' for two groups which satisfy the conditions of the 
proposition. We have 

GAr_i(iV, 1, -1) = G^^N-i){N, 1, 1), with 6(i) = -J (mod N). 
G2{N,n,-l) = G-2{N,n,l), with e(j) = -j (mod iV = 2n+l). 
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Gn-i{N,1,~1) = G2{N,n,-l) for N ^ 2n + l. This is the isomorphism 
between the synimetries of the niaximal chains (see 15.3.11) and the symmetries 
of the Eight with an odd number of bodies. Here ki] = —1, krj' — — n, p = 
n - 1, 6(j) = (2n - l)j = -2j (mod N). 

6*3(4, 1, —1) = Gi(4, 1, 1). This is the isomorphism between the symmetries 
of the Gerver solution (see l5.3p and the symmetries of the relative equilibrium of 
the square. Here krj = —1, hf ~ 1, p = 1, 6{j) = 3j (mod 4) = —j (mod 4). 

63(5, 2, 1) = Gi(5, 1, 1). This is the isomorphism between the symmetries 
of the 3 lobes chain for 5 bodies (see ....) and the symmetries of the relative 
equilibrium of the pentagon. Here krj — 2, ki]' — I, p — 1, 6(j) = — 3j 
(mod 5) = 2j (mod 5). 

Remark The difference between conditions 2) and 2') is most easily under- 
stood on the example of the groups Gi(4, 1, 1) and Gi(4, 2, ±1). The horizontal 
relative equilibrium of the square possesses both symmetries but the two ac- 
tions are not isomorphic: indeed, relative equilibria in an inclined plane are 
symmetric under Gi(4, 1, 1) but not under the symmetry group of the Hip-Hop 
Gi(4, 2,±1). This comes from the fact that condition 2) is trivially satisfied 
because p = while condition 2') is not. 

4 Local continuation 



After reduction by translations and rotations, the (67V— 10)-dimensional reduced 
phase space is the sum of a (4iV — 6)-dimensional "horizontal" subspace and a 
(27V — 4)-dimensional "vertical" subspace, both invariant under the linearized 
flow. The eigenvalues of this flow are that of the variational equation, the only 
difference being that the generalized eigenspace corresponding to ±iLUi has lost 
four dimensions, two horizontally (corresponding to fixing the value of the an- 
gular momentum and to rotations around it), and two vertically (corresponding 
to rotations around a horizontal axis). 

We have already studied the vertical part of the spectrum. A thourough 
study of the horizontal part appears in [Moj . See figures [5] and [B] Note that 
purely imaginary horizontal eigenvalues exist only for A^ > 6 bodies. 
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Figure 5: Horizontal spectrum for three to six bodies 



32 




* * ,** + ** + **** 



Figure 6: Imaginary parts of purely imaginary horizontal eigenvalues of the 
regular iV-gon, as functions of A^ e {7, ..., 30} are upper bounded by aj„ (solid 
line) . 



Recall that for each vertical frequency ojk the vertical eigenspace, most of the 
time, has dimension 4. The only exception is when A^ = 2n is even and k = n; 
if moreover no other vertical or horizontal frequency is an integer multiple of 
Wn, Lyapunov's classical theorem implies the local existence and uniqueness of 
the Lyapunov family. In the other cases, one can try to get existence from the 
following theorem of Weinstein, here in the form given by Moser, of whom we 
have kept the notations. 



Theorem 19 ( [Mos', p. 743]) Let H be a Hamiltonian on i?^" with an equi- 
librium point at the origin, and C be the linearization of the Hamiltonian vector 
field at the origin. Assume that R^" ^ E ® F , where E and F are invariant 
subspaces of C such that all solutions of C in E have the same period T > 
while no nontrivial solution in F has this period. Moreover, assume that the 
Hessian D^H(0) restricted to E is positive definite. Then, for sufficiently small 
e, on each energy surface H{z) = H{Q) + e^ the number of periodic orbits of H 
is at least ^dimE. 

We will see in the next section that the theorem applies to the vertical 
frequencies provided that no unexpected resonance occurs. The cases N = 3 
and N = 2n with Hip-Hop symmetry are respectively studied in |CF2j and 
|BCPS| (see also [MSj for the case of the regular n-gon with a central mass at 
its center). 

4.1 Partial convexity of the energy 

In order to apply the Weinstein-Moser theorem to some vertical frequency tui, 
one must check that the energy level sets, restricted to the space of 2TT/ujg- 
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periodic solutions of the linearized vector field, are compact in the neighborhood 
of the relative equilibrium. 

Let Vi be the vertical eigenspace of the frequency we {1 < I < N/2), and 
Til be the plane tangent to the homographic motions. (Recall that the total 
eigenspace of wi contains Vi ©Hi.) If no other frequency, horizontal or vertical, 
is an integer multiple of (possibly equal to) uJi, in order to apply the Weinstein- 
Moser theorem it suffices to prove that the energy is convex on Vi ® Hi for £ = 1 
and on V^ for 2 < £ < N/2. Below we will prove the fact that the quadratic 
part H is positive definite on the whole vector space 



^ = Hi ( 



h<e<N/2> 



The first cases, studied in the following section, all satisfy the required 
non resonance condition. Furthermore, numerical experiment suggests that the 
purely imaginary horizontal eigenvalues, in general, cannot resonate with the 
vertical eigenvalue iw„, at least, for they are smaller in module (see figure [6l 
where the horizontal frequencies have been computed using the factorization of 
the stability polynomial which is described in |Mo| )rl When this is true, the 
proposition below and the Weinstein-Moser theorem show in a weak sense the 
local existence of Lyapunov families associated with cli„, in particular. 

Proposition 20 After reduction by rotations, the restriction of the quadratic 
part of the energy to T is definite positive. 

We first prove two lemmas. The frequency Ui plays a special role, for its 
eigenspace always includes a horizontal plane, namely the plane tangent to the 
homographic motions. 

Lemma 21 After reduction by rotations, the restriction of the energy to Vi®Hi 
is definite positive. 

This lemma allows to apply the Weinstein-Moser theorem [Mosi p. 743] to 
the frequency uji provided that the reduced linearized equation have no solution 
of period 2-K /lo\ outside Vi 0Hi. It is a straightforward generalization of |CF2i 
lemma 2.1]. 

Proof. For the same reason as in jCF2[ lemma 2.1], the lift of \>\ ® Hi to the 
non-reduced phase space is tangent to the submanifold 

£i = {{x, x) e (M^)^ X (R^)^, 3p,(j>0 3R, S e SO3, x = pRC, x = aSK} 

where, as before, C is the central configuration (1,^, ..., C^~^) and J is the 
horizontal-rotation operator. It suffices to prove that the restriction of H to £1 
is definite. For the sake of completeness, we shortly explain the computation. 



■^Note however that the largest horizontal frequency seems asymptotic to uin, when N tends 
to infinity. 
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Before reduction, {p,a,R,S) e (M+)2 x SO{3) x SO{S) are (generalized) 
coordinates on £i and the restriction of H is 



H\.. = T E iMc^ii'-7 E 



1 



0<i<JV-l '^ 0<j<fe<Ar-l "^ ^ " 

As in 11.31 the mass dot product by C of Newton's equation shows that 

1 



Hence 



0<j<k<N~l "^ ^ ' 



We compute the reduced system by first quotienting by the fuU group 50(3) 
and then fixing the length of the angular momentum L. This amounts to re- 
placing {p, a, R, S) by (p, a, R~^S) and imposing the relation 

N-l 



pa||5](C^0)Ai?-l5(^c.lC^0)||:=||L| 



i=0 

Any element of a neighborhood of the Identity in 5*0(3) can be uniquely written 
as exp A, where A is an antisymmetric 3x3 matrix. In particular, 

-c b 
R-^S = exp I c -a 
-b a 

(b^ + c^) ah ac 

ah ~{c^ + a'^) be 

ac be — (a^ + b^), 

where the dots represent terms of order higher than or equal to 3 in a, b, c. 
Locally, the angular momentum equals 

N Mo/ ^ W + b^ + 2c^ 





\L\\ = pCTLUlN 



l--{a^ + b^+ 4c2) + 



In the coordinates (a, b, c, d = a — 1) on the space £i reduced by the rotations, 
the reduced Hamiltonian is 

^{-l + a^ + b' + 4c' + d^) + ..., 
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which proves the lemma. 

We will now study the other eigenspaces Vg. Usually Ve is 4-dimcnsional, 
generated by the sum of the Lyapunov "cylinders" S{N,i,±l). As we wiU see 
in examples, the first interesting case is that oi N ^ 5, £ = 2. The case where 
A^ = 2n is even and ^ = n is special: Ve = V„ is then 2-dimensional (and the 
Lyapunov center theorem applies with no need of convexity property, provided 
that no other frequency is an integer multiple of Cij„). 

Lemma 22 If 2 < £ < N/2, Vg lies in a submanifold of fixed vertical angular 
momentum and in the reduced space the restriction of the energy to Vi is definite 
positive. 

Proof. A point on Vg can be uniquely represented by (x, x) with 

Xj = {cosj6,smj9,acosj£9 + bsmj£9) 

Xj = {—uiismj9,ujiCOSJ9,cu!isiTij£9 + dujiCosj£9), 

with 9 ~ 2tt/N, as before, and where a, b,c,d G M; (a, 6, c, d) is a coordinate 
system on Vg, except in the case N = 2n and £ — n, where (a, c) alone is a 
coordinate system (with b = d = 0). 

The statement on the angular momentum comes from the fact that the 
angular momentum of the j-th body is 

( cos j9 \ / —uji sin j9 

sin j9 \ '^ \ ^1 cosj^* 

a cos j£9 + b sin j£9j \ujfc cos j£9 + LOgd sin ]££ j 

its vertical component depends only on the radius of the horizontal circle, and its 
two horizontal components boil down to degree-one trigonometric polynomials 
in j{£ + 1)9 and ]{£ — 1)9, whose sums with respect to j e {0, ..., A^ — 1} are 
zero if £ > 1 . 

We now need to compute the quadratic part of H with respect to a, 6, c, d. 
A straightforward computation shows that the kinetic energy is 

hence positive definite with respect to c and d. On the other hand, the potential 
part depends only on a and b and, at the second order. 



-U 



E 



1 



0<j<k<N-l "^'^ ^■'' 

>r- 1 f 1-a^ sin^j + k)f - 6^ cos2(j + k)f + 

" ^ 4| sin(i - fc)f I V 2a6sin(j + fc)f cos(j + /c)f + • • • 



„ .r-^ (a sin(j + k)^ ~b cos(j + k)^ 
= -Null + 2_] -i ^ ^2 v-' ^2 



4|sin(j-fc)f| 
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which is positive, definite if ^ ^ N/2 and null \i (. — N/2. 

We now complete the proof of the proposition. 
Proof. We want to show that in the neighborhood of the Lagrange relative 
equilibrium the quadratic part of the energy H is positive definite on the space 

T = 1-Li ^ ®l<t<N/2^l- 

The two lemmas above show that the quadratic part of the energy H is 
positive definite in restriction to Tii ® Vi and to each V^ , t>2. As eigenspaces 
of a symplectic linear operator (namely, the linearization of the reduced vector 
field) are orthogonal, it follows that the aforementioned spaces are pairwise 
symplectically orthogonal. Since additionally they are invariant by the flow of 
the quadratic part Q of H\jr^ Q does not have off-diagonal terms. Hence it is 
positive definite. 

4.2 The first cases 

We will now describe the simplest families. It is interesting that new qualitative 
features appear for each successive value of N from 3 to 6. 

A natural name for each Lyapunov family we study would be the label 
S'(iV, fc, rj) of the tangent Lyapunov cylinder in the inertial frame. On the other 
hand, for the sake of concreteness we like to look at a given family in rotat- 
ing frames where some fixed symmetry Gr/s(-/V, fc,??) holds. This amounts to 
considering the Lyapunov family as a deformation of the Lyapunov cylinder 
Sr/s{N, k, rf). Among the infinitely many possible choices, we have selected the 
symmetry of the most remarkable periodic orbit of the (global) family in the 
inertial frame. Of course, this is a subjective choice and a family may bear 
several names. 

Three bodies The case of A^ = 3 bodies is fully analyzed in |CF2| . which we 
now summarize. The eigenvalues of the linearized reduced vector field (which 
can be computed following |Mo| ) are: 

Horizontally: zbiwi, (±-\/2 ± «) (-Ji 

Vertically: ±jwi 

The eigenvalue ztzwi has multiplicity 2 (1 horizontal and 1 vertical), corre- 
sponding to two Lyapunov families: 

- The horizontal family of homographic motions; as a curiosity, see figure [7] 
for a view of a member of the family in a choreographic rotating frame. The 
so-called proper degeneracy of the Newtonian potential shows up here in that 
the family has no torsion. 

- The vertical Pi2-family, tangent to the Lyapunov cylinder 5(3, 1, —1). 

The Pi2-faniily (S'(3,l,— 1)) Because of the 1 — 1 resonance between 
the vertical and horizontal eigenvalues iiwi, one has to understand the third 
order normal form of the vector field in order to prove that there exists a unique 
Lyapunov family tangent to a non-trivial solution of the (WE); in an adequately 
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Figure 7: Homographic solution in a choreographic frame, rotating twice per 
period 



rotating frame, this unique family possesses a G2(3, 1,— 1) = -De symmetry, 
which is the Fi-symmetry described in [CFMj . See figure [H As elsewhere, the 
solutions have been computed using 

- the Dopri-8 algorithm [HNWj, as implemented in the C language by F. 
Joutel and M. Gastineau (Astronomic et Systemes dynamiques, IMCCE), for 
integrating Newton's equations 

- a hybrid Newton algorithm, implemented in the Gnu Scientific Library 
|GSL| . for solving symmetry equations. 

The germ of the Pi2-family was described in [Malj as the family of rela- 
tively periodic solutions of the three-body problem with the highest possible 
symmetry; it is immediately after he learnt about the existence of the Eight 
that Christian Marchal noticed |Ma2| that these solutions turned into chore- 
ographies in a rotating frame. This is the origin of |CFlj . where we give its full 
scope to this remark, and prove in particular that: 

- the flow of the reduced three-body problem in the neighborhood of the 
Lagragange relative equilibrium has a unique 4-dimensional center manifold A/", 

- the energy levels within J\f are 3-spheres, 

- the restriction of the flow to an energy level has a Poincare section diffeo- 
morphic to an annulus, whose return map is a monotone twist map which is the 
identity on the boundary corresponding to the homographic motions. 
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Four bodies 



Horizontal eigenvalues 

('±iv/56-14\/2± i\ wi ~ (±0.8595325038 ± i) ui 



(±y42V-16 + 18V2-49 -(±0.6394812009 



± i Y42 \/-16 + 18 V2 + 49) ^ ±0.9533814590 i)wi 

Vertical eigenvalues 
±iuji =±§v/2V2 + l 
±iw2 = ±iv%wi ~ ±1.2155625241 iwi 

There are now two non-trivial vertical Lyapunov families. 

The 4-body, 3-lobe chain family (S'(4, 1,-1)) The Weinstein-Moser 
theorem applied to the reduced vector field and to the eigenvalue ±10^1 (using 
lemma [21I and the fact that uji/uj2 — 2^^!^ ^ N) implies the existence of at least 
one more Lyapunov family corresponding to this frequency in addition to the 
horizontal homographic family. 

Numerical computations show that the Lyapunov cylinder 5(4, 1, —1) is tan- 
gent to a family which is choreographic in a rotating frame which starts making 
two full turns per period in the negative direction (see the top part of figure [T5|). 
This does not follow from the Lyapunov theorem, for the frequency Lo-y is also 
the frequency of the horizontal homographic family. 

Yet without uniqueness it is theoretically not obvious that the family shares 
the 62(4, 1, — l)-symmetry with iS'(4, 1, —1), in the rotating frame making one 
turn per period, at the limit of the regular pentagon. The proof of the local 
existence and uniqueness of this family would require similar computations (of 
a normal form of the third order) and arguments as for the Pi2-fainily described 
above. 

The 4-body Hip-Hop family (S'(4, 2, ±1)) Since ll>2 has multiplicity one 
and since the only other frequency, wi, is not an integer multiple of (^27 it follows 
from Lyapunov's theorem that there exists a unique family of relatively periodic 
orbits bifurcating from the square with vertical frequency close to uj2- The 
eigenmode iS'i(4, 2,±1) is tangent to the family and, by uniqueness, the family 
shares its symmetry. Hence the family is the Hip-Hop family, with symmetry 
group Gi(4, 2, 1) = Z/2Z x Z/4Z (see figure H). It is studied in [BCPSj . 
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Five bodies 



Horizontal eigenvalues 

~ (±0.9391304549 ±i)wi 
~ (±0.8281366700 ± 0.8822431635 i)u;i 
~ (±0.5028535236 ± 0.8822431635 ijwi 
Vertical eigenvalues 



±ZC^l = ±^ ^ 



V2V5-Vf(V5+2) 
1/4 \/5-3/4+2 (l/4 VE+l/^y 



±ZW2 = ±iWl'\/2 < 



V^\/5-V5(ll/4+3/4 Vs) 
V2V5-V5(V5+2) 



±1.3281310261 iwi 



The 5-body, 4-lobe chain family (S'(5,l,— 1)) Figure [TT] depicts it in 
rotating frames which start by making three full turns in the negative direction 
(symmetry 64(5, 1, —1)). 

Contrary to its analogue with four bodies, the maximal vertical frequency 
072 now has multiplicity two. Using lemma [22] and the above expression of 
UJ1/LU2 ^ N, the Weinstein-Moser theorem proves the existence of at least two 
Lyapunov families. Numerically, there are exactly two families, which thus 
display the same symmetries as their corresponding Lyapunov cylinders: 

The 5-body Eight family {S{5, 2, -1)) It owes its name to the 6*2(5, 2, -1)- 
symmetry it acquires in a frame with r/s = 2 (figure fT2|l . 

The 5-body, 3-lobe chain family (S'(5,2,l)) It owes its name to the 
non-maximal chain 6*3(5, 2, l)-symmetry in a frame with r/s = 3 (see figure [H]) . 

Six bodies 



Horizontal eigenvalues 



Vertical eigenvalues 



±1102 
±iW3 



±m/^±^ 



V3 
3 



±i\J3 + 
±fyi7 



V3 



±0.4687282051 iwi 
±0.4211614102 wi 

(±0.9499800584 ± 0.8151022048 z) wi 
(±0.2986755303 ± 0.8151022048 z) wi 
(±0.9893611078 ± i) lui 



±1.3991678967 iwi 
±1.5250481798 icji 
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From six bodies on, some horizontal eigenvalues are purely imaginary, which 
could give rise to horizontal Lyapunov families. For example, could the chore- 
ography shown in figure 3 (4th of first column) of [5] belong to such a family? 
Moreover the purely imaginary horizontal eigenvalues can resonate with the ver- 
tical frequencies, on which we focus (see the remark before lemma \n\ however ] . 

The 6-body, 5-loop chain {S{6, 1,-1)) It is the 6-body instance of max- 
imal chains. Its proof of existence has the same status as with 4 and 5 bodies. 

Some 6-body Hip-Hops Since 6 is not prime, the 6-body problem has a 
richer set of Hip-Hop families than problems with fewer bodies. 

The Weinstcin-Moscr theorem and numerical experiments indicate the exis- 
tence of a Hip-Hop family sharing the Hip-Hop symmetry of and tangent to the 
Lyapunov cylinder 5(6, 1, —1); it belongs to an invariant problem having as few 
dimensions as the two-body problem, where bodies are symmetric with respect 
to each other within two groups of three (see part 5 of figure [4]). 

The second frequency gives rise to two Lyapunov families S{6, 2, 1) and 
S'(6,2,— 1), each belonging to an invariant problem having the dimensions of 
the three-body problem, where bodies are symmetric with respect to each other 
within three groups of two (see part 6 of figure 2] for 51(6, 2, —1)). 

Since all other frequencies (vertical or horizontal) are smaller than lu^, by 
the Lyapunov theorem the third vertical frequency gives rise to a unique Hip- 
Hop family sharing the symmetries of S'(6,3,±l). In a frame rotating once 
per period (at the limit of the hexagon) , this family is both a Hip-Hop (bodies 
are symmetric with respect to each other within two groups of three) and a 
partial choreography (bodies chase each other by pairs, along three distinct 
closed curves). 



5 Global continuation 

We are interested in the following two questions: 

- Existence: Does the range of frequency rotation w of the frame over which 
the family exists contain 0? 

- Uniqueness: Can one take the frequency n? as a monotonous continuous 
parameter over the whole family, i.e. has the torsion constant sign? 

Minimization under the G^/s (./V, fc,7/)-symmetry is a natural tool for the 
existence question. It will turn out that, because of isomorphisms between the 
actions of different groups Gr/s{N,k,rj) (cf. section [375|l . its use is essentially 
restricted to the case k = n, i.e. to the largest vertical frequency. 

As a preliminary study we apply the results of section 12.21 to estimate inter- 
vals of the rotation frequency of the frame over which the relative equilibrium 
family itself is the sole absolute minimizer under the Gj./s{N, fc,77)-symmetry. 
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5.1 Minimization properties of the A^-gon family under 
Gr/s{N, k, r7)-symmetry 

Consider solutions x{t) ~ e^^^y{t) of Newton's equations which, in a frame 
rotating with frequency tu, become s-periodic loops y{t) € A*^, where G = 
Gr/siN,k,ri). A relative equilibrium solution with the required symmetry ac- 
complishes one turn in time - in the rotating frame. Hence it has frequency 
LOi = 2TT- + 137 in the inertial frame. 

Referring to proposition [8l we now study conditions under which the inf of 
the positive A's for which there is a solution x{t) — e^^*y{t) with y{t) £ A'~^ of 
the equation 

—X = XAx 

is equal to 1, where A is computed with the mutual distances fij of the relative 
equilibrium of frequency 2t:- + zu in the inertial frame. 

Let ^(i) be defined by x{t) — ^(\/At). As ^(t) is a solution of — ^ = A^, each 
component of ^(t) is of the form J^k '^feS""'^'''*? where uji,uj2, ■ ■ ■ , ojk, ■ ■ ■ are the 
frequencies of the (WE) and 

iOk = (2.TT h Tjj). 

UOX S 

The linearity of the differential equation and the structure of the action of 
Gr/s{N, k, rj) allows us to study separately solutions lying in the horizontal plane 
and solutions lying on the vertical axis: if y(t) is invariant, so are the loops h(t) 
and v{t) of A^-body configurations in R'^ obtained from y{t) by projecting each 
body on the horizontal plane and the vertical axis respectively. 

1) VERTICAL SOLUTIONS. A solution x{t) of the (WE) of the form 

(xo(i),--- ,XN-i{t)) withVj e Z/NZ,x,{t) = (0,0,4(0), 
can be written 

x'^{t) = ^Re UuiC' + viC')e'^"^'] , j E Z/NZ. 
I 

The symmetry conditions impose not only s-periodicity but even 1-perio- 
dicity. Indeed, if we choose f] = 0,S = 0,q ~ 1, we get 6 = I E R/sZ, hence 

a:|(i)=y|(t)=2/|(t-l)=x^^(t-l), 

that is 

^Re {uiC^^+viC^){l-e-'^"^)e"^"^' =0. 
I 

It follows that for any /, {mC^'^ + viC^'-){l - e"*'^'^) = 0, that is 
either ui = vi = or 3m e Z, lji\^ ~ 2TTm. 
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The second possibility occurs for a single index I provided no relation of the form 
maji' = m'uji exists with m, m! € Z. In any case, as we are only interested in the 
values that A can take, we can restrict the attention to solutions corresponding 
to a single index I such that uJi^/X = 27rm. 

Now, the invariance of x{t) under the action of Gr {N, k, rj) means that for every 
{9, 6, P, C) € Z/sZ X Z/2Z X Z/2Z x {-1, +1}, such that 6* - § - krjj^ = g € Z, 



Re 



(uzC-" +i';C-")e*"'^* - (u(C^^-'+*^' ^y^^iU+^nyii-^l^+'^iV^iit-e)] 



= 0. 



Choosing S — 0, (3 = 1,^= 1, hence 9 = 1/2, wc get that to = 1 + 2p is odd. 
Choosing 5 — 0, /3 ~ 0,£_ ~ —1, hence 9 = 0, we get that for all j E Z/NZ, 
{ui — ui)C,^^ + {vi — vi)C,-'^ = 0, which implies that ui and vi must be real. 
Finally, for ^ ~ 1 the condition is that for every j, S, 



uiC^ (l - e^^< ^^) j + viC'^ [1 - e'^H T^)j 



= 0, 



which implies that one of the two coefficients ui and vi must vanish and that, 
according to which one does vanish. 



I = ±mkri mod A^ = ±(1 + 2p)kri mod N. 
Finally, the inf of the corresponding VA's is the inf of the 

(1 + 2p)27r 



(l+2p)27i- 



that is 



inf ■ 



Ull 



P ^(l+2p)fer, 



2tt'- 



where the frequencies W(i_|_2p)feT; = ^-{i+2p)krj involved are the ones of bifurcating 
vertical solutions with the required symmetry (in particular, ujkri — ^k is always 
among them). One deduces immediately that the condition inf A > 1 reads 



vj + 27r- 



< y = inf 



Ul 



P>0 UJ(i+2p}krt 



|(l + 2p)2^|. 



where the retriction to the p > comes from the fact that if p' 
has |1 + 2p'\ = |1 + 2p\, hence uJi+2p'kr, = tJi+2pfer,- 



-{1+p), one 



2) HORIZONTAL SOLUTIONS. The general "horizontal" solution 

c-(eo,---,eA'-i), withe,(t) = (e"w,o). 



of — ^ = A^ is such that, identifying M^ with 
such that 



there exists kQ,ai,bi,ci,di G C 



Vj, i-it) ^ko + J2 [(°'e^"'* + bie-''^'')C' + {cie''"' + die-''"')^ 



E 

1^0 



If A^ = 2n + 1 is odd, to each of the n pairs A; = A^v-; = —(^f of non-zero 
eigenvalues of the matrix which defines A, corresponds an eight dimensional 
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space of solutions parametrized by the four complex constants ai,bi,ci,di. This 
gives the dimension An of the phase space after quoticnting by the translations. 



If A^ = 2n is even, to the eigenvalue A„ corresponds only a four dimensional space 
of solutions because (" — ("^ ^ —1. In this case the two complex parameters 
are ai + q and bi + di. 

We shall now study the conditions imposed on horizontal solutions by the in- 
variancc under Gr (N, k,ri). 

1) The condition of s-periodicity on y{t) imposes that, in addition to a)o = 0, 
only one frequency uji is present: indeed, for all j and all i, 

= ^[(e*"'" - l)(a,C^' + QC^'')e*"'* + {e"^" - 1)(6jC^'' + diC'^)e"^''], 

where 

ai ~ ujivX — n7, /?; = — a);v A — tu. 

Hence, as uji ± lji/ ^ ii I ^ I' , 

Vj, I, (e^"'^ - l)(a,C^"' + qC^"') = (e^^^^ - l)(6zC^' + diC') = 0, 

As aiC^^ + ciC^'- (resp. kC^^ + diC,'^) cannot be equal to for aU j = 0, • • • , iV - 1 
except if a; = c; = (resp. bi ~ di =0), one deduces that there exists a 
(necessarily unique) I ^ such that e*"'* = 1 or e*'''* = 1, that is 



ais 



— {(liivX — zu)s ~ in2TT, or /3(S = (— a)(V A — zu)s = m2TT, 



I.e. 



yj{t) = iko + {aiC'+ciC^')e'^^'^\0) or y,{t) = (fco + (fei^' + djC''')e*'^^*, 0). 

Both cases are identical hence we suppose that the first identity holds. Hence 
we are looking for the invariance of 



y,{t)^{ko + iuiC^'+viC')e'^-'^',0), 



corresponding to 



VA: 






vj + I-k"^ 



n7 + 27rf 
This means that, for all (6^, b, (3, £,) e R/sZ x Z/NZ x Z/2Z x {-1, +1} with 

9 kri — = Q G Z, 

2 'n ^ ' 



^27ria^,^ 



one has for aU j e Z/NZ and for all t, yj{t) = e'^'^'°'f^^._^_g-^{^{t - 6)) 
a ~ ^6 — -^ (mod 1). As ^^ = 1, this is equivalent to the following: 



where 



kn + {uiC'' + viC'')e 



'jl\ i2ir^t _ 



fco* + (flR(-''+')' + ««C^^+^)')e2-^(*- 



^27via 
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Let us first take ■f = 1 and hence uj — ui and vj — vi. The condition becomes 
the vanishing of the expression 



„27r»(f+a-^; 



#+^)('l_e2"(-#+"-T^) 



27ri(- 



for all j and all t. This implies the vanishing of one or the other of the following 
complex numbers : 



1 



„27ri(-^;+Q_ 



or B = 1 



^2mi-^l+a- 



9) 



- First case: ^ = 0, that is yj{t) = (fco + u/C-"e2"'?*,0). 

- Second case: B = 0, that is yj{t) = (fco + viC^'e'^'"-'?\0). 

In both cases, fco = unless a = mod 27r but fco does not play any role in 
the determination of A. 

Taking now ^ = —1, we see that ui in the first case, vi in the second one, 



must be real. Replacing a by 
S 



4f mod 1, we get 



T — - Tfl n T — Tfi 

,,(/-!) + 0eZ or -^il + l) + ^eZ. 

JS s N s 

Finally, replacing 6* by ^ + krjjT + g, we obtain the following conditions: 



l-l 



2 ' ""IN 

r — m 



-krj 



5 r — m ,13 , „ 



or 



^{l + l) + kr] 



S r — m ,0 , „ 



Taking P ^ 0,q = Q,S = 1, this implies that 

I _ I _l_ r^zinji^ = mod N in the first case, 

— (/ + 1) + ^^^^^krj = mod N in the second one. 

Taking now /? = 1, this implies in both cases that, for any q G Z/sZ, 

r — 771 , 1 



-( — h 9) = mod 1. i.e. 3p G Z, m — r — 2ps. 
s 2 

It follows that 1 = 1 — 2pkri mod N in the first case, I = —(1 — 2pkri) 
mod N in the second one. Notice that both cases correspond to the same type 
of solutions 

y,{t) = (fco + a,C^"(i-2pMe2"(7-2p)t^o), 

(and that, when in addition 1 — 2pfc77 is prime to N, these solutions are up to 
translation and scaling, a relative equilibrium of the regular N-gon with a differ- 
ent period and a relabelling of the bodies, that is one of the relative equilibrium 
solutions found in section 1X5)) . This shows that the case where A = B = Q will 
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not lead to different solutions and hence need not be studied separately. In both 
cases, as toi = uj-i, 



Vx = 



UJi 



^l-2pkri 



vj + 27r-j — 4p7r 



w + 2t:'- 



Notice that the value p — corresponds to A = 1. Hence, in order to insure 
that the minimum value of A is 1, it is enough to insure that A is always > 1. 
This amounts to the following conditions, for which one should remember that 
the values of p must be such that 1 — 2pkri ^ mod N . 

i)li p> Q and w + 27r- > 4p7r, the condition is 

r 
i^i-2pkrt - wi)(tu + 27r-) < -4p7rwi, 

which can never be satisfied because for every I, cui — uji > 0; this implies 

that 

r 

w + 2tt- < Att; 

s 



a) If p > and < tu + 27r- < 4p7r, the condition is 

r 

{uji + uji^2pk-n)(.^ + 27r-) < ipnuJi; 

in) If p > and ru + 2tt- < 0, the condition is 

r 

{uji - uji^2pkrj){-^ + 27r-) < ipnuJi; 

iv) If p < and zu + 2tt- < Api:, the condition is 

r 

{uji - uji^2pk-n)(.^ + 27r-) < Apiruji, 

which can never be satisfied for the same reason as above; and this implies 

that 

r 

w + 27r- > -47r; 
s 

v)li p <0 and 4p7r < ro + 2tt- < 0, the condition is 

r 

-(wi + tJl-2pfer,)(ra7 + 27r-) < -ApiTLJi; 

vi) If p < and w + 27r- > 0, the condition is 

r 

(.^i-2pk-n - uji){taJ + 27r-) < -Apiruji. 

Finally, the horizontal conditions are that —H— < w + 2tt- < +-ff+, where 
if,=inf/47r,inf*^ ^^""^ ^^""^ 



P>0 \CJi + ijJi-2pkri ^l+2pkr] — ^1 
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if_=infU7r,inf*' ^^^'"^ ^^'^"^ 



p>0 ywi + U)l+2pkri ^l-2pkr] — ^1 

where the inf * means that the index k of any ujk involved can never be 0. 
Recalling the vertical condition 

{l + 2p)27r\, 





m + 2Tr- 
s 


<y = 


= mf 

p>0 t^(i+2p)fe, 


we obtain the 









Theorem 23 The following condition implies that the relative equilibrium solu- 
tion of the equal mass regular N-gon with frequency 27r^ in a frame rotating with 
frequency zu is the sole absolute minimizer of the action among paths which in 
the rotating frame are s-periodic loops with the Gr/s{N,k,ri)- symmetry of some 
solution of the vertical variational equation: 

- inf(F, H-) <w + 2-K- < inf(F, H+). 
s 

When k — n, the horizontal contribution disappears; in this case, the corre- 
sponding Lyapunov families can be searched for as absolute minimizers of the 
action with the given symmetry constraints: 

Corollary 24 The following condition implies that the relative equilibrium so- 
lution of the equal mass regular N-gon with frequency 27r- in a frame rotating 
with frequency tu is the sole absolute minimizer of the action among paths which 
in the rotating frame are s-periodic loops with the Gr/s{Nin,vi)- symmetry: 

-V <vj + 2tt- < V. 
s 

Proof. As LOi < oJn, it is enough to prove that 

'^ >!- and '' > ' 



^l±2pn ± Wi UJn Wl=F2pn ± ^1 ^n 

in all the cases where no ujq is implied. But all these inequalities are implied by 

2uJn > 0Jl±2pn + ^1- 

In the following sections, where the numerical side plays an important role, 
we first study the two families associated with the highest frequency w„: the 
Eight families when N is odd and the Hip-Hop families when N is even. The 
first ones generahze Marchal's Pi2-family ([CMl lMa2l [CTMl [CFTl |S]) which 
links the Lagrange equilateral relative equilibrium to the Eight, the second ones 
generalize the family which links the relative equilibrium of the square to the 
original Hip- Hop ([CVl|TV]). These are precisely the families for which global 
minimization of the action under the Gr/siN, k, 77)-symmetry constraint may be 
used for global continuation. But even then, no unicity of minimizers is proved; 
hence the continuity of the families so obtained is in question, even if global 
topological continuation results can be used. 
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We then study the famihes associated with the lowest frequency uji, which 
lead to Simo's chains with the maximal number of lobes (A'' — 1 lobes for A^ 
bodies): only numerical results are given here because these families are only 
local minimizers of the action. This is related to the fact that, starting with 
TV = 4, there exist isomorphisms between the actions of some of the groups 
Gr/s{N,k,ri) (see section [X5|) . We also give examples of chains with a non 
maximal number of lobes, where similar phenomena do occur. 



Remark It follows from corollary[T7]that absolute choreographies are dense in 
the Lyapunov families as soon as torsion is present (i.e. as soon as tu does vary 
along the family). This results from the fact that a path wich is Gr/s{N, k, rf)- 
symmetric in a frame rotating with frequency vj is Gr/s+ro (-^j fcj^)-symmetric 
in the inertial frame. (In the action of the group, this corresponds to the re- 
placement of a by a + vjO.) 

5.2 Families associated with the highest frequency u;„ 

They satisfy V — 2tt^^^ and it follows that the relative equilibrium remains the 
unique global minimizcr as long as 

-27r— !- <w + 2tt < 2ti—-, 
i.e. until one reaches the vertical bifurcation of the corresponding family. 
5.2.1 The Eight families: G = G2(2n + 1, n, -1) 





Figure 8: The P12 family. Filled and hollow circles represent bodies at times 
i = and t = T/12 



When n = 1 one gets the original Pi2-family of Christian Marchal (see fig- 
ure[5]). In this case, the minimization property of the relative equilibrium family 
was already proved in |BTj under the weaker assumption of the choreography 
symmetry. Note that for planar solutions, 27r is replaced by 47r. We leave to 
the reader the pleasure to check that this was a priori obvious. 

When n > 2, aproof of the existence of a G2(2n-|-l,n, —1) = £'2n-symmetric 
Eight is still missing. The existence of a Z3„-symmetric Eight is proved in (FTj 
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but, while highly probable, the fact that it is automatically D2n-symmetric is 
not proved. Figures fT2l and [T8l depict the case n = 2. 



5.2.2 The Hip-Hop families: G = d (2n, n, ±1) 

The original Hip- Hop (figures P and [TB]) corresponds to the case n = 2. For all 
values of 71, the existence of the Hip- Hop family is proved in [TV] , with the usual 
proviso that no uniqueness, hence no continuity, is proved. It seems very likely 
that the natural end of the family is a pair of simultaneous n-tuple collisions 
for w = ui2, but this not proved either. The heuristic explanation is that the 
phenomenon is the same as for a fixed center pb with angle a = 27r: in the 
inertial frame, each body must turn by exactly 27r during time T/2). Moreover, 
all along the family, the minimizing solutions should possess the brake symmetry. 




Figure 9: The original Hip-Hop family (Hip-Hop in the middle) 



Remark Illustrating the remark at the end of 15. 1[ figure [TC] shows the two 
simplest spatial choreographies (in the inertial frame) in the Hip-Hop family for 
n = 2. First described by S. Terracini and A. Venturelli, they correspond to the 
symmetry groups G3(4, 2, ±1) and G5(4, 2,±1). On the other hand, figure [TO] 
shows the original Hip-Hop in a frame rotating with its frequency W2- 





Figure 10: The Hip-Hop in a frame rotating with its very frequency 
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5.3 Chains families: the role of angular momentum 

The chains are planar choreographies of the equal-mass iV-body problem, similar 
to the Eight but with a number £ > 2 oi lobes. The first example with £ > 2 is 
Gerver's chain with 3 lobes and 4 bodies [CGMS| . Because of their symmetries, 
chains with an even number £ of lobes have a vanishing angular momentum 
while this is not the case when £ is odd. This fundamental difference accounts 
for a totally different behaviour of the corresponding Lyapunov families. We 
have studied numerically the cases of 4 and 5 bodies and have observed: 

- when £ is even, a complete unfolding of the corresponding unchained poly- 
gon into a vertical chain; 

- when £ is odd, a complete unfolding into a horizontal chain (in a frame 
which is still rotating) , followed by a plane family which continues to a horizontal 
chain. 

Due to the existence of isomorphisms of the actions of different symmetry 
groups Gr/sii^i k, rj), global minimization of the action may succeed only for the 
Eights (£ = 2). 

5.3.1 Maximal chain families: G = Gn^i{N, 1, -1), A^ = 2n -I- 1 

When observed in a frame which rotates A^ — 1 times per period in the negative 
direction, the non-trivial solutions of the (WE) corresponding to the frequency 
uji give rise, to an infinitesimal choreography which unfolds the A^ — 1 circles. 
When TV = 2n -f 1 is odd, the full family hopefully continues up to the fixed 
frame into a vertical zero angular momentum planar chain with A^ — 1 lobes. 
We recall that it was proved in section [33] that the action of Gn-i{N, 1, —1) is, 
for A^ = 2n-|- 1, isomorphic to the one of the Eight with A^ bodies G2{N, n, —1). 
We study now the simplest case where the maximal chain family and the 
Eight family differ. 

5.3.2 The 4-lobe chain and the Eight for 5 bodies: two isomorphic 
actions of the symmetry groups 

Four-lobe chain: G ~ 64(6, 1, —1) When A^ = 5, we have checked numeri- 
cally that, indeed, the 6*4(5, 1 — 1) family continues up to a vertical planar chain 
with four lobes. During the unfolding, the two central lobes become smaller and 
flatten more rapidly while the the two exterior lobes remain during a long time 
almost vertical before flattening down as the Big Ears of Big Brother. 
We have 

T/ = 27r, H+^in — — — , iJ_ = inf f 27r, 47r — — — 1 == 27r. 



Hence, on the left hand side {m + 8tt < 0), the estimate goes all the way till 
the bifurcation of the chain family but on the right hand side, it stops before: 
—2tt<w + 87t< Att ^ "^^ . This is due to the existence of the Eight family (see 
figure fT2)) which, as we have recalled, has the same symmetries up to reordering 
(see figure [HI) . 
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Figure 1 1 : The family of the 5-body chain with 4 loops 

Five-body Eight: G = G2(5, 2, -1) The solutions of the (WE) corespond- 
ing to the 2 families with frequency W2 (recall that they differ only by the 
replacement of C,^ by C,^) give rise to chains with respectively 3 lobes and 2 
lobes, the order of the bodies in the choreography being respectively the one of 
the retrograde and the direct stellated pentagon. We are interested now in the 
second one. 
We have 



F = 27r 



UJ2'' 



H. 



inf (47r, Stt- 



Ui 



UJ2 — LUl 



127r- 



LUi 



UJl +UJ2 



) = 47r, i?_ = 8tt 



UJl 



UJl +UJ2 



It follows that A^ = 1 as long as — 27r— < ru + An < 27r^^, i.e. until we 
reach the vertical bifurcation of the Eight family (see figure [T5)) . 




Figure 12: The family of the 5-body Eight. 



Remarks. 1) The case of G = G4(5, 1,— 1) is the first one where _ff+ < V. 
This is because the second frequency uj2, while ruled out from the G-invariant 
vertical solutions, is present into the G-invariant horizontal solutions. 

2) The estimate H+ yields a lower bound of the action of the member of 
the Lyapunov family which bifurcates at n7 = (—2^2 — wi)^^. In particular, 
at this point, the action is higher than it is at the bifurcation point, because, 

< 



UJ2 



2^1 



Ui+ijJ2 ' 



5.3.3 The 3-lobe chains for 4 or 5 bodies 

Four bodies When A^ = 4, the non-trivial solutions of the (WE) correspond- 
ing to the frequency uji, more precisely to the group G3(4, 1,±1), give rise to 
3 lobes chains when one starts with a frame which rotates two full turns in 
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the negative direction. When the rotation decreases, the chain starts opening 
but the central lobe decreases and when it reaches approximately —0.9 turns 
by period, an almost collision occurs. After that, the solution flattens to the 
horizontal plane (with vj not yet vanishing) where a second bifurcation leads to 
the Gerver solution (figures [T3] and [T7|) . 




Gerver' s choreography 
Figure 13: The planar Gerver family as a secondary bifurcation 



Five bodies An analogous scenario is observed for N = 5 for the frequency 
UJ2, more precisely with the group 6*3(5, 2, 1). but in this case, the family never 
gets close to collision (figures [14] and [T9| and the central lobe is bigger than the 
two extreme ones. 




Absolute choreography 

Figure 14: The 5-body 3-lobe chain family and its planar secondary bifurcation 
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5.4 Action diagrams 

Figures [T5Hroi display the action of families of solutions as a function of the 
frequency w of the rotation of the frame. Each figure corresponds to a fixed 
symmetry in the rotating frame (recall that the isomorphisms between group 
actions are associated with relabelling of the bodies as described in thcorcm [^5|) : 

Figure[l5] G2(3,l,-1) = G_2(3,l,l) 

Figureim G'i(4,2,±l) 

Figureini G3(4,l,-1) = Gi(4,l,l) 

Figurelia G4(5,l,-1) = G2(5,2,-l) 

Figureim G3(5,2,l) = Gi(5,l,l) 
The solutions are represented in perspective in R^. This holds in particular 
for solutions lying in the horizontal plane (relative equilibria and horizontal 
secondary families in figures [17] and [T9|) or in a vertical plane (Eights or chains 
in figures [15] and [18]) . 

Display of the implications of theorem 1231 Fat segments on the ro-axis 
and the corresponding fat part on the graph of the action indicate intervals 
on which the relative equilibrium of the regular 7V-gon is the unique absolute 
minimizer for the given symmetry in the rotating frame. It appears that the 
theorem detects global phenomena i.e., the presence of different branches of 
Lyapunov families. 



Symmetry with respect to n7 = of figure llSI This symmetry corresponds 
to the isomorphism G2(3, 1,— 1) = G-2(3,l,l). Analogous symmetries could 
have been shown on the other figures. 

On the w < 0-half (resp. w > 0-half) of the figure, the configuration of 
the relative equilibrium is a positively (resp. negatively) oriented equilateral 
triangle. In the rotating frame, it makes two turns in the positive (resp. nega- 
tive) direction. Notice that in the incrtial frame the four possible combinations 
of orientation of the configuration and direction of motion occur at the four 
bifurcation points — Swi, —uji, toi and 3uji. 

Also, in figure [TS] the continuity of the Lyapunov family implies a change of 
sign of the a;-coordinate of the body when m goes through 0. This is consistent 
with the Gr{N, k, 77)-symmetry: if a loop is Gr{N, k, 77)-symmetric, the same is 
true for its image under the rotation of angle tt around the vertical axis. 

Cases without torsion The GAr_i(iV, 1, — l)-Lyapunov family starting at 
w = —Nuji (resp. +Nloi) from the relative equilibrium which completes minus 
(resp. plus) one turn per period in the inertial frame has no torsion: indeed, 
it is merely the family obtained by rotating the horizontal relative equilibrium 
around the y-axis. As this family has constant action, this fact accounts for the 
end of the action diagrams in figures [T^] [13 and [TS] 



About figure 1181 The family of solutions which in figure [18] bifurcates at 
vj = {—2uj2 — oJi)^^ is, up to scaling, symmetric of that which, in figure [TOI 
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bifurcates ai zu — toi ~ Scji. More precisely, the latter family transforms into 
the former one by: 

1) reflection about zu — 0, which corresponds to the isomorphism 6*3(5, 2, 1) = 
G_3(5,2,-l); 

2) translation of —5lu2 along the ro-axis, which transforms the G_3(5, 2, — 1)- 
symmetry into the symmetry group G2(5, 2, —1) of figure \TE\ (see the remark at 
the end oflSTTj): 

3) scaling by lui/uj2- 

In the rotating frame, the family starts as the Pi2-family but, at some point, 
two more loops develop on the supporting curve which then flattens to a hori- 
zontal planar (relative) choreography with two loops, described twice per period. 
The figure below shows this (quasiperiodic) solution in the inertial frame. 




54 




Figure 15: The action of the P12 family and of two times Lagrange solution in 
the rotating frame 
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Figure 16: The action of the original Hip-Hop family (including two absolute 
choreographies discovered in |TV| ) 
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Figure 17: Action of the 4-body, 3-loop chain family, and of the planar Gerver 
family 
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Figure 19: Action of the 5-body, 3- loop chain family, and of the corresponding 
planar family (cf. [S]) 
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6 Appendix: Fourier expansions and the torsion 

Approximate Fourier expansions of Gr/s(-/V, fc, 77)-symnietric Lyapunov families, 
in the same spirit as Marchal's computations in [Mali IMa2] , allow to evaluate 
the torsion under a regularity hypothesis. 

6.1 The symmetry ansatz 

We are looking for local one-parameter families of solutions of the iV-body 
problem which, in a family of frames rotating with frequency 

r 

vj — uji (jJk + n7, 

s 

are periodic of period T = s—. We will suppose that ujk — 2tt, and hence T = s. 
Such solutions are of the form 

x,{t) - {h,{t),z,{t)), h,{t) = e^(-i-^"'^+^)*^,(i), J = 0,- • • ,7V- 1, 
with 

h,{t)= Y^ afe^^.it^ z,(t) = Re(^ fofe^^'^^*)- 

l — ~oo 1— — 00 

Moreover, we ask the solutions in the rotating frame Xj{t) = {hj{t), Zj{t)) to be 
symmetric under the action of Gr/si^, k,r]) described in section 3. Recall that 
an s-periodic loop of configurations x{t) — [xi{t), • • • , XN{t)) is invariant under 
the action of Gr/s{N, k, rf) if and only if, for every (6*, (5, /3, ^) G G2 representing 
an element of Gr/s{N, k, 77), i.e. such that 9 — ^ — kij-^ = m G Z, one has 

yj e z/NZ, x,{t) = px^(j+s) {i{t - 0)), 

where the action of p on M^ = C x M is defined by 

p(h, z) = (e'^'^^/i*, e"'^^) with a = -e (mod 1), 

s N 

where /i? = /i if ^ = +1 and h^ ^h\i^^ -1. 

Translated in terms of Fourier coefficients, this reads: 



. ^ ^.2.(a-ie) £0-+5) 



for the horizontal components and 



Re {h^ - e"P%^'+'^ ) = and bj = e*(-/5?-2-^e);,«0-+^) if / ^ 
for the vertical components, implying restrictions on the coefficients aj and bj : 
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Horizontal coefficients We have 

/ r-l 5 r-lj , fr-l \ S 

a- -9= 6 - — = (^ + m + kri - 1 — ' 

s s N s ^2 ' \ s ' J N 

Fixing ^ and S and changing /3, the angle e*^'^^"^"^) takes at least two different 
values as soon as ^-^ is not even. Hence 



r 



-I 



af = if is not an even integer. 

s 

If i-di = 2p is even, the symmetry conditions become 



Applying these identities with both values ^ = 1 and ^ = — 1 , we get by difference 



VPj Ji 'j'r — 2ps — '^1 — 2psi nence Vp, a^— 2ps — ^r~2ps' 

r-2ps ^^6 '"sal and 

V„ o „■?' - -J27r(2pfer,-l)i 



So, the coefficients a°_2ps ^^^ '"^al and 



Finally, as the center of mass is at the origin, we have that for all /, X]i=o '^l ~ ^^ 
hence 

'Af-l 
^^-.2.(2pfc,-l)^|„0_^^^^0^ 

which implies 

al_2pg = for all p such that 2pfcr7 —1 = (mod N). 

To summarize, 
a^ — unless possibly for I = r — 2ps with 2pkri — 1 7^ (mod N) 

V„ o „■?' - -J27r(2pfcr,-l)i 

Vertical coefficients Giving successively its two possible values and 1 to 
/?, we get that 

Vj, Re6^o = 0. 

Moreover, 

— £ 9 = — (£ ) — m kri 

takes two different values for /? = and /? = 1 as long as £ — - is not an even 
integer, which implies that, for I ^ 0, bj = can be different from only if I is 
of the form / = (2p + l)s, with p an integer. 
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Now, a I ^ 0, choosing 6 = — j and giving its two possible values ±1 to <^ 
and noticing that (3^ = (3 (mod 2), we get 

hence 6° G M. Finally, as the center of mass is at the origin, V^, X]i=o ^i = ^^ 
that is 

(n-1 \ 

gi27r(4(l-i)-mi) y^ g»27rfc7)i^ 6? = 0. 

Finally, bj* (and hence all 6;^) must vanish unless k-ql ^ (mod sN). 
To summarize, 

bj* = unless possibly for I = {2q + l)s with (2g + 1)A;7? 7^ (mod N), 

V9,J, "(2q+l)s - ^ "(29+1)^ 

In particular, 

^^z^*) = X! "^29+1 cos 27r(2g + l)(t + fc77— ), 

where we have used the notation 

C29+I = b{2q+l)s + ^-(2<,+ l)s' 9 = 0, 1, • • • 

For example, these conditions for the P12 family (that is for the group 
62(3, 1, —1)) become 

«[> e M, af == if either I = 1 (mod 2) or ^ == (mod 3), 

50 e M, 6/ = if either I = (mod 2) or ^ = (mod 3), 
which coincide with the conditions found by Marchal in [Mal| . 

6.2 The regularity ansatz 

In the rotating frame, the solutions wc are interested in are the ones tangent to 
the cylinder of solutions of (VVE) 

(AoC^e'2'^5*, eReC''''^e*2^*), j = 0, • • • , TV - 1, 

where as usual ( = e*^ . 

In the inertial frame, these solutions are of the form 
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hj{t) 



^i{uji+^)t\^ n 



P 



pi27r[-2pt-(2pfcr,-l)i] 



2jW =Xl"=2«+iCOs27r(2g + l)(t + fc?7^) 



g>0 

where the coefficients a^(e) and ci{e) are real and the summations are respec- 
tively over all integers p G Z and q G Z+ such that 

2pkri - 1 7^ (mod N) and (2q + l)k-q ^ (mod N). 

We will assume that, in the reduced phase space where the relative equilib- 
rium becomes an equilibrium, the one-parameter family of s-periodic solutions 
we are looking at generates an analytic foliation of an analytic 2-dimensional sur- 
face, more precisely, that it is the image under an analytic embedding of the triv- 
ial planar family ii(i) — ecos27ri/s, v(t) — esin27ri/s. As (cos''' 27ri/s)(sin ^-Ktjs) 
has no harmonics of order larger than k + I, this implies in particular that the 
Fourier expansion of the coefficient of e" in x{t) does not contain harmonics of 
order larger than n. This justifies taking 

as our main local parameter and making the following ansatz: 



a. 




r—2ps 



^2pe 



|2p| 



0(e 



|2p|+2N 



C2,+i = C2,+ie^'+^ + 0(e^«+^) 



,29+3 \ 



which in the case of the P12 family coincides with the one made in [Mal| . 

Using both ansatz, an identification in the equations of motion allows us to 
determine the leading coefficients a'l and cj, as well as the dominant term in w, 
in terms of e. The dominant terms in the components of Zj(t) are 

Zj{t) = ecos27r(t -|- krjj-) + Cse^cosGnit + krjj-) + O(e^), 



while those of hj (t) are 



E 

0=0, ±1 



,0 , 

1 — 2ps 



a. 



'^<-^pt+W)+Oie^), 



that is 



hi = e'('^i+''" )* 



j^ i2^(-2t+q' 



[{ 



^ t2^{2t+i 



+ Oie^), 



where we have used the following notations (we have directly set to zero the 
coefficient of e in a^ and w because this is an immediate consequence of the 
equations in the next section): 
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The Eight families G'2(A^ = 2n + 1, n, -1, ) 



The Hip-Hop families Gi(iV = 2n,n,±l) 

h,{t) = e*(-i+-)* (a°e*2.^ + a0e*2.(2*+(2"+i)^) + „o ^g«2.(-2t-(2„-i)^) ^ . . .^ ^ 

The maximal chain families G'2n(2n +1,1,1) 

/^j(i) = e^(-i+-)* (a0„e*2.^ ^ ^0^^^g«2.(2t-^) ^ ^0^_^g,2.(-2t+3^) ^ . . . j ^ 

except if A^ = 3 where the last term is absent. 

6.3 Identification of dominant coefficients 

It is remarkable how well the two ansatze above make it possible to identify 
Fourier expansions of the Lyapunov families: the regularity ansatz makes the 
system of equations block-triangular, while the symmetry ansatz makes the 
solution unique for a given choice of group. 

Notations 

1 j — I 

Aji = A^pju Bji = — sin2 27r(fc77^^), 6^,;^ = O^i - 0j[pj,[p] . 

Finally, we recall the restrictions on the values of p and q which come into the 
sums; in particular, X]p=±i ^^^^ mean the sum restricted to those values of 
p — ±1 such that 2pkri — 1 7^ (mod N). 

Mutual distances With the notations above, 



-ii-Kpt- 



p 



^0 -^i^Yl '^29+1 

Mod 0{e^) we obtain 



7 / 

cos27r(2(j+ l){t + k'qj-) - cos 2TT{2q + l)(t + krj—) 



\\hj-hi\\^^{AoPjif + 2Ao 



P=±i 
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z, -zi 



j I 

cos27r(t + ki]—) — cos27r(i + kr]—) 



e\ 



that is 

\\hj~hi\\^ = {AoPjif+2AQPji 



a + ^ A2pPj[p]i[p] cos47r(p< + 0ji - 6'j[p]j[p]) 



6^ 



|z^- - z,|2 = 4sin^ 27r(t + k-q-^^^) sin^ 27r(fc77:^ )e2, 



which, using the notations above, becomes 



\h,-hi\\' = Al 



1 + ^— a + X! ^2pPj[p];[p] cos47r(pt + Q^i-p) e^ 
^J' V p=±i / 



\zj ^zi\^^ 2A,iBji ( 1 - cos4^(t + fcr;^) ) e\ 



2N 



Finally, ||xj — xjH '^ equals 






1- 



3 / a + Ep=±i-42pPj[p]([p]Cos47r(pi + ej7;p)+ \ ^ 
^i \ B,i (l - cos47r(t + fcr;|^:^ '' 



0(6^ 



We will now plug in the obtained expressions into the equations of motion 

hi - hj .. v^ zi- z, 



^. = E- 



It^J 



\Xl Xj 



13' ^J 






Horizontal equations Recall that 



\Xl Xj 



|3- 



p=0,±l 



hence 



. -s --^ / ^0 „i27r(-2pt+%l)„ 



0{e^ 



that is 



j[p] ^ 



- Y. A^p [-(^1 + 76^)2 + 87rp(^i + 76^) - IGttV] e*2.(-2p*+^),2 + ^(,4)^ 
p=±i 



j[p] ' 



Ep=±i ^2p [-c^? + Sttp^i - 167rV] e*2.(-2pt+2fi) 
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Finally, as 



=±1 



the horizontal part of the equations splits into the 0-th order equation, which 
is nothing but the equation satisfied by the relative equilibrium: 



i^, (^o^J-') 



and the second order equation 



-c.i(2Ao7 + c.ia)e^2.^ ^ ^ ^^^ ^_^2 ^ g^^^^ _ ^^^2^2] g.2.(-2pt+2W) ^ 



p=±i 



;^j \ p=±i 



2pe ^^''^'n^jM 



^4 / a + Ep=±i^2pPj[p];[picos47r(pt + ej/;p)+ \ ^ 

We will call {Hj) this equation. 

Now, because of the symmetry ansatz, the N equations {Hj) are equivalent to 

one of them, for example (ifo)i which is of the following form 

U + T/e*'"'* + We-''^''' = 0, 

where the complex numbers U, V, W are affine functions of the 4 real unknowns 

a, A2, A_2, 7- 

Moreover, it turns out that, because of the invariance of the expressions Aqi, Bqi, ■ 
under the change of I into —I, the coefficients of U,V,W are indeed real. As 
(Ho) has to be satisfied for all values of t, it is equivalent to the three real affine 
equations 

U = V = W = 0. 

Vertical equations At the order of approximation 0(e*), we have 



hence 



zi - z. 



Zj{t) = e(cos27r(i + fc?7-^) + C3e''cos67r(i + /cr7-^) + 0(e^). 



Uji sin27r(t + krj'-——) + Caw^/e^ sin67r(t + krj——) + ©(e'^) 
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where 

Uji = ~2am2TT{kr]^^), Vji = -2 sin 671(^77-^^). 

The vertical equation (V^) is obtained by identifying the e^-terms in the identity 
(mod 0{e^)) 

-47r2 cos 27r(i + krjj-) - SGn^Cse'^ cos67r(t + krjj-) = 
^\\xi-Xj\\'^^ (ujis[n2'K{t + kT]——) + CsVjie'^ sine-Kit + kj]-——) + 0(e'*) 

The identification of the terms of order give the vertical variational equation 
(VVE) which is already satisfied by Zj{t) = cos27r(i + k-qj^). 
The identification of the terms of order 2 gives the two remaining relations 
beween the five unknowns a, A2, ^-2, C^, 7. They are of the form 
Re (Xje^^TTt-) ^ Q^ j^g (y,e*6'^*) = 0, 

where Xj is an affine function of a,A2,A_2 and Y,- is an affine function of 
^21^-21^^3- As above, the symmetry ansatz implies that the N equations 
(Vj) are equivalent to one of them, for example (Vb) and the invariance of the 
expressions Aqi, Bqi, ■ ■ ■ under the change of I into —I implies that Xq and Yq 
are real, hence that (Vq) reduces to exactly two real equations. 

Identification of coefficients Finally, we have 5 real equations which are 
affine in the 5 unknowns 

a,A2,A^2,C3,j. 

They can be solved in the following order (assuming non degeneracies which 
will prove true in the cases investigated below): 

- In (Hq) the term in e" depends only on Aq, which it determines. 

- In (Hq) the term U in e^ constant with respect to time depends only on a 
and 7 and can be used to eliminate a. 

- In (Hq) the coefficients V and W of £2g±i4irt q^^q.^^ ^q eliminate A±2. 
(For the Eight families, it can help to remember or it can be checked, that 
a° = A2e^ + ... = 0.) 

- In (Vb): the coefficient of e^gieTrt (^f^gj- simplification by e) allows to elim- 
inate C3. (In the case of the P12 family, C3 = 0.) 

- In (Vb), the coefficient of e^e^'^'^* eventually allows to compute 7. 

Remark The fact that, in the rotating frame, the horizontal period is half 
the vertical one follows from the invariance of the problem under the symmetry 
with respect to the horizontal plane (materialized in the action of P). This 
fact accounts for the generic aspect displayed by the solutions in a frame which 
accompanies the rotation of the regular A^-gon, as exemplified on figure 1101 
Notice that the same type of behaviour, with the same explanation, is observed 
in the restricted problem for the vertical solutions originating from the Lagrange 
points (see |Zaj). 
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6.4 Torsion of the first cases 

Recall that we are computing the torsion of the Lyapunov families which bifur- 
cate at w — ~-LiJk + ^1- Since wi > 0, 7 is necessarily > 0. 

Three bodies 

The family of rotated Lagrange solutions (5*2(3, 1, 1)) 

7 = as expected. 

The Pi2-family (52(3,1,-1)) 7 = — (6^^)'^^ ~ 16.589 

This value is consistent with Marchal's computation [Malj . Indeed, writing 

with primes quantities introduced in his book, he finds that the frequency shift 

of the Pi2-family is 

a' , 3 ,9 , 

- + h.o.L = -cf^h.o.t.^ 

to be compared with our 76^. But cijj = 1 while uji = 2it. Hence 

2-^cf = 7^^ 

(that Marchal's masses are not the same as ours does not come into play in this 
equality between frequencies). Besides, his c'-^ is the amplitude of oscillations of 
zi — Z2 at the first order and the edges of the equilateral triangle he considers 
have length 1. On the other hand, our e is the amplitude of oscillations of zi 
and we consider a triangle whose edges have length aoA/3, with a^^ = 4V37r^. 
Hence 

whence the above expression of 7. 
Four bodies 

The 4-body, 3-lobe chain family (^2(4, 1,-1)) 

^- 41x2W/3 (^^-^^)(^^+0 >9 
The 4-body Hip-Hop family (S'i(4, 2,±1)) 
4692 7r^/3 
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2446096 X 2^4^2^2 + 1 
Five bodies 
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440V2 + 981) > 19 



The 5-body, 4-lobe chain family (5*4(5, 1, —1)) 

6 2i (17995 - 7823 V5) tts 

7 = ^^ -1 > 5. 

5981 X 53 

The 5-body Eight family (5*2(5,2,-1)) 

_ 7ri(l + V5)i/g(l + 3V5)i/g 
^~ 211 X51/3 X 1155/6 761809^ 

(230356305630646272 V5 + 954953942246092800) VTH^ V375^~T ^ > 19 
+10427244637440119808 V5- 3532108049365632000 J 

The 5-body, 3-lobe chain family (53(5,2,1)) 

/5 + l)^(3V5-l)* 



69129994912 x 5i/3 x lin/e 

/ (522541527 V5- 177004875)- \ 

i (11543868 v^ + 47855700) V\/5 - 1 V3 \/5 - 1 J ^ 



Six bodies 



The 6-body, 5-loop chain (55(6,1, 1)) 

48 TT^/"^ 

7= —(32634^/3- 39889) (15-4^3)2/3 >3 

' 781199x69621/3^ '^ ' 

A 6-body Hip-Hop (5i(6, 2, -1)) 7 > 14 (long expression) 
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